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We develop a systematic analytic approach to the problem of branching and annihilating random 
walks, equivalent to the diffusion-limited reaction processes 2 A — > and A — > (m + 1) A, where 
m > 1. Starting from the master equation, a field-theoretic representation of the problem is derived, 
and fluctuation effects are taken into account via diagrammatic and renormalization group methods. 
For d > 2, the mean-field rate equation, which predicts an active phase as soon as the branching 
process is switched on, applies qualitatively for both even and odd m, but the behavior in lower 
dimensions is shown to be quite different for these two cases. For even m, and d near 2, the active 
phase still appears immediately, but with non-trivial crossover exponents which we compute in an 
expansion in e = 2 — d, and with logarithmic corrections in d — 2. However, there exists a second 
critical dimension d' c ps 4/3 below which a non-trivial inactive phase emerges, with asymptotic 
behavior characteristic of the pure annihilation process. This is confirmed by an exact calculation in 
d — 1 . The subsequent transition to the active phase, which represents a new non-trivial dynamic 
universality class, is then investigated within a truncated loop expansion, which appears to give a 
correct qualitative picture. The model with m = 2 is also generalized to N species of particles, 
which provides yet another universality class and which is exactly solvable in the limit N — > oo. For 
odd m, we show that the fluctuations of the annihilation process are strong enough to create a non- 
trivial inactive phase for all d < 2. In this case, the transition to the active phase is in the directed 
percolation universality class. Finally, we study the modification when the annihilation reaction 
is 3 A — > 0. When m = (mod 3) the system is always in its active phase, but with logarithmic 
crossover corrections for d = 1, while the other cases should exhibit an directed percolation transition 
out of a fluctuation-driven inactive phase. 
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I. INTRODUCTION 

Reaction-diffusion processes, in addition to being of direct physical interest, present a relatively simple class of non- 
equilibrium stochastic processes in which the role of fluctuations, so important to the study of equilibrium statistical 
physics, may be investigated Q|. Although such processes have been studied for many years through various kinds 
of mean-field or self-consistent approximations, more recently there has been some progress in taking into account 
the effects of fluctuations through the use either of exact methods or of systematic renormalization group theories || , 
backed up by extensive numerical simulations. 

It has become clear that, once the constraints of detailed balance are abandoned, these systems may exhibit a 
greater richness of phenomena than their counterparts close to equilibrium, leading, for example, to new types of 
critical behavior characterized by new universality classes. It might appear that the freedom allowed by giving up 
detailed balance would lead to such an unacceptably large number of possibilities as to make the concept of universality 
of limited usefulness, and indeed, from the point of view of the proliferation of microscopic cellular automaton models 
this would seem to be true. However, as with equilibrium phenomena, the introduction of a coarse-grained description 
through a field-theoretic formalism makes it possible, at least in principle, to make a universal classification feasible. 

All such coarse-grained descriptions of stochastic processes, whether or not they satisfy detailed balance, divide the 
dynamics into two parts: a deterministic evolution corresponding to the mean-field dynamics, and the noisy effect of 
fluctuations ignored in such a description. When detailed balance is present, the form of the noise is prescribed, so 
that the field-theoretic version of the coarse-grained description is narrowly determined. In the absence of detailed 
balance this is not so, yet it may be shown that the particular form of the noise chosen is crucial to determining the 
universal behavior near a non-equilibrium critical point. One way to avoid this ambiguity, which works very well for 
reaction-diffusion processes, is to develop a field-theoretic description of the problem which is exact, being directly 
derived from the underlying master equation, and to perform coarse-graining, in terms of a derivative expansion, only 
at a later stage. This method, developed some time ago by Doi and by Peliti (see also Ref. Q), relics on the use of 
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a Fock space formalism very similar to that employed for quantum many-body systems. More recently, field-theoretic 
renormalization group methods have been systematically applied to the dynamic held theories describing some simple 
processes to obtain what is, in principle, a complete characterization of their universal properties |6| ||]. 

It is the purpose of this paper to describe such a field-theoretic investigation, largely but not entirely based on 
renormalization group methods, into a class of reaction-diffusion processes which are characterized by their description 
in terms of branching and annihilating random walks. (A very brief account of this work has been presented in Ref. ||.) 
These processes were originally studied in the 1980s by Grassberger et al. jlO| as exceptions to the apparent rule that 
systems exhibiting non-equilibrium transitions from a trivial absorbing state into a non-trivial noisy 'active' state 
should be in the same universality class as directed percolation (DP). While this terminology originally refers to a 
realization of this universality class describing percolation processes which are asymmetric with respect to a privileged 
space dimension |liq ], more commonly it now applies to a dynamic universality class in which this dimension is 
interpreted as time. In terms of reaction-diffusion processes, an example of such a system would be one with a single 
species of particle A, undergoing diffusive behavior, single particle annihilation A — > 0, and branching A — > 2A. There 
is always a trivial absorbing state, with no particles and consequently no fluctuations. For sufficiently low branching 
rate, this is the only stationary state, but for larger values of this rate, another, non-trivial 'active' stationary state 
appears. The DP transition occurs at the point where this state first appears. The robustness of the DP universality 
class has been tested by measuring its critical exponents in numerical simulations on many different types of system, 
and may be understood in terms of the simplicity and irreducibility of the dynamic field theory which describes it 
|l3|| . Indeed, this 'Reggeon' field theory was first studied in a completely different context of particle physics. Among 
the non-equilibrium phase transitions described by this universality class are those occurring in the contact process 
[ 14 1 . the dimer poisoning problem in the Ziff-Gulari-Barshad model [l5[ , and in certain autocatalytic reaction models 

The model studied by Grassberger et al. p0| ] was a rather complicated probabilistic cellular automaton in one 
space and one time dimension, with antiferromagnetic interactions in which the 'particles' actually appear as kinks or 
domain walls. It was observed that, while this model appears to exhibit a transition from an inactive, absorbing state 
to a non-trivial active state, the various critical exponents are far from those of directed percolation. The authors of 
Ref. Jl^] already realized that this may be a result of the fact that these kinks are conserved modulo 2 in the basic 
dynamics, but a subsequent field-theoretic analysis of this problem appeared to shed no light on this discrepancy: 
rather it suggested that the system should, once again, be in the DP universality class, indicating that there was 
apparently no way of taking into account the constraint of local 'parity' conservation in a field-theoretic treatment 
p| . One of the main points of our analysis is to show where this argument went wrong, and why such systems, at 
least in low dimensions, should be in a new universality class. 

Subsequent work showed that similar effects occur in other systems with such a conservation law, and suggested 
that, at least in one dimension, they all fall into this new universality class, characterized by branching and annihilating 
walks with an even number of offspring |l7] 2(J. These are reaction-diffusion systems with the underlying reaction 
processes 2A — > and A —>■ (m + l)A, with m even. In one dimension, when m = 2, another realization of this 
universality class occurs in certain classes of dynamic Ising models which violate detailed balance [^0| , pr|j2ll , where 
once again the 'particles' correspond to domain walls. In this case, the transition from the inactive phase (no domain 
walls in the asymptotic stationary state) to the active phase corresponds to the disordering transition of the Ising 
spins [ p2|p3| . In addition, other realizations of this universality class have been found and also been investigated 
numerically [§4]]2f|. 

On the other hand, when m is odd and there is no conservation law, one expects to recover the DP exponents, and 
this appears to be the case, at least in one dimension fl0| , |l7|j2"6| , p7f . In higher dimensions, the numerical situation 
is less clear. Takayasu and Tretyakov, in rather short simulations, found no evidence for a non-trivial transition for 
d > 2, for either m even or odd p7| . This is certainly consistent with the picture obtained from mean-field theory, 
as we shall see. But the situation in two dimensions is, we believe, more subtle. It is just this dimension at which 
the fluctuation effects become important in the pure annihilation process 2 A — > 0, giving logarithmic corrections to 
scaling, and they cannot therefore be ignored there. In fact we find in our analysis that these fluctuation effects are 
sufficient to drive the existence of an inactive phase for small but finite values of the branching rate, at least for 
odd m, and hence the existence of a subsequent transition in the active phase, which is in the two-dimensional DP 
universality class. 

For m even, on the other hand, the fluctuations are not effective until a second, lower critical dimension which 
we estimate to be « |, and which is certainly above d = 1. The existence of two critical dimensionalities make a 
systematic e-expansion type of approach to this problem infeasible. However, we have performed a truncated loop 
expansion which has the correct fixed point structure to describe the new d = 1 universality class. One striking feature 
of this is the existence of dangerous irrelevant variables which make the identification of the various critical exponents 
in terms of RG eigenvalues slightly delicate. In fact, we may only make this identification strictly in the active phase, 
and our analysis leaves open the possibility of violations of the scaling relations between these exponents and those at 
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the critical point and in the inactive phase, which have been written down on the basis of simple scaling arguments. 

Our aim, therefore, is to investigate the following competing reaction processes of identical particles A, which 
perform a random walk with diffusion constant D: 

annihilation : k A — > , rate A& , (1-1) 
branching : A — » (m + 1) A , rate a m , (1-2) 

with k,m integers, k > 2, and m > 0; '0' denotes an inert state. Occasionally, we shall in addition allow for a 
spontaneous creation of particles in pairs, 

pair production : — > 2 A , rate r . (1-3) 

The corresponding mean-field rate equation for the average density n(t), neglecting any local fluctuation effects, 
reads 

dn(t) 

—^ = -k\ k n(t) k +ma m n(t) + 2T, (1.4) 



with an initial density n(0) = no. If we set the pair production rate, which acts as an 'external field', to zero, Eq. (1.4) 
has the stationary solutions n = ('inactive state') and 

/ma 

T = ° : n = "- S (w) (L5) 



('active state'). Provided the 'mass' ma m is positive (and no > 0), the solution of Eq. (1.4) approaches the asymptotic 
density n s exponentially for t — > oo, 

n 

r = 0: n(t) = - - rT/Jk^T) ' ( L6 ) 

f 1 + [(n s /no) fe_1 - 1] e-^- 1 )™^*] 

On the other hand, if o~ m = 0, this exponential long-time behavior is replaced by a power-law decay, 

u m = t = : n(t) = U ° -TJTuZT) ■ ( L7 ) 

[1 + fc(fc - 1) A fc rig" 1 1] /( J 

In this sense, <j TO = cr c = may be viewed as a critical point, described by the following set of critical exponents 

r = 0, a m = a c : n{t) cx t~ a , (1.8) 
r = , a m I c c : n s oc (cr m - cr c ) /3 , (1.9) 

characterizing the long-time decay at the critical point, and the growth of the asymptotic density in the active phase, 
respectively; according to Eqs. ( p..5| ) and (1.7) their mean-field values are ao = /3q = l/(fc — 1). Taking into account 
the diffusive propagation of local density fluctuations, two more independent exponents may be introduced which 
describe the divergence of the characteristic length and time scales as the critical point is approached, 

T = 0, a m — » a c : £ oc \a m - o- c \~" , (l-!0) 
t c ot t z * \a rn - a c \- z » , (1.11) 
a m = a c ,T^0: £ oc |r|-" T , t c oc |r|-*^ , (1.12) 

with the mean-field values vq = 1/2, v t q = l/(d + 2), and zo = 2. 

These mean-field exponents of course originate in the scaling dimensions of the parameters D, a m , and r of 
the model. Introducing a momentum scale k, and measuring times in terms of lengths squared, 

[f] = «T 2 , [n] = n d , (1.13) 



(as the particle density scales as an inverse volume in d dimensions), Eq. (1.4) implies the 'naive' dimensions 

[D] = k° , [r] = n d+2 , (1.14) 
[<r m ] = n 1 , [A fe ] = K 2 -( fe -D d . (1.15) 
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Therefore, both particle production processes appear to be relevant perturbations to the pure annihilation reaction, 
and for a m = r = the annihilation rate becomes dimensionless at the (upper) critical dimension d c — 2/(fc — 1). 
Thus, one expects the mean-field rate equation (|0|) to provide an at least qualitatively correct description for d > 2 
in the case k — 2, and for d > 1 in the case k = 3. Because of the relevant perturbation a m , the situation at the 
critical dimension already requires a careful analysis, and certainly for d < d c fluctuation effects ne ed t o be taken into 



account seriously. E.g., for the pure annihilation process with k = 2 it is known exactly that Eq. (1.7) is replaced by 
the slower asymptotic decay 

k = 2, d < 2 : n(t) oc r d/2 , (1.16) 
due to strong particle anticorrelations which cannot be smoothened out by the diffusion in low dimensions [pp|§] ■ At 



d = 2 for k = 2 and d = 1 for k = 3, one finds logarithmic corrections to Eq. (1.7). However, for k > 4 the above 
rate equation should be essentially correct in all physical dimensions. In the bulk of this paper, we shall therefore 
be concerned with the case k = 2, and only briefly discuss the three-particle annihilation reaction (k = 3) combined 



with branching in Sec. VIII 



The outline of the rest of this paper is as follows. After introducing the field-theoretic description of these processes, 
with an emphasis on the subtleties which occur when correctly taking into account processes in which a conservation 



law is satisfied, in Sec. Ill we look close to the critical dimension d — 2 of the pure annihilation reactions to discover 
whether the branching processes are relevant (thus leading immediately to the active phase) or otherwise. For m 
even we find that they are in fact irrelevant close to two dimensions, leading to a transition into the active phase at 
zero branching rate, with, however, non-mean-field exponents which may be evaluated within an e expansion. This 
calculation suggests, however, that this situation might be reversed in sufficiently low dimensions, and in Sec. |y| 
we confirm this by an exact calculation in one dimension, using both field-theoretic and exact lattice methods. The 
subsequent non-trivial transition into the active phase is described within the truncated loop expansion in Sec. 
and the various subtleties connected with the unusual fixed point which emer ges are investigated. Comments are 



given on the feasibility of extending this calculation to higher orders. In Sec. VI a generalization to N species is 
introduced and solved, in the hope of elucidating the case N = 1. However, it turns out that the behavior as soon as 
N > 2 is quite different, leading to yet another universality class. In the next Section we turn to the case of m odd, 
and use diagram resummation and RG methods to argue that for d < 2 the fluctuation effects drive the exis tence 



of a non trivial inactive phase, counter to the predictions of mean-field theory. After a brief discussion in Sec. VIII 
of the corresponding set of problems when branching is added to the three-particle annihilation process 3A — » 0, we 
conclude with a summary of our results and possible directions for future theoretical and numerical work in this area. 



II. MASTER EQUATION AND FIELD THEORY 



In o rder to include fluctuation effects systematically in a mathematical description of the diffusion-limited reactions 
(1.1)— (1.3), we first write down the corresponding master equation and then derive a field-theoretic representation for 
it, following standard procedures ||-||,|l(J . The latter may then be treated with standard methods like perturbation 
theory, and the long-time and long-distance scaling behavior may be inferred via the application of the renormalization 
group. 

Master equations for the temporal evolution of the probability distribution P(a; t) for a configuration a have the 
general form 



dP(a;t) 
dt 



E^- 



(2.1) 



where R a ^p denotes the tra nsit io n rat e for the process a — ► (3. For random walkers on a lattice (with lattice constant 
&o) subject to the reactions ([iTl|)-(|L"3|), any configuration is characterized by the integer site occupation numbers n,, 
i = 1, . . . , N, with m — iVyi(t), and the master equation may be decomposed into the balance equations for the 
diffusion, annihilation, branching, and pair production processes at each site i, 



dP{{n r };t) 
dt 



E 



dP(m;t) 



dt 



dP(m;t) 



D 



dt 



dP(m;t) 



dt 



dP(m;t) 



dt 



(2.2) 



with 



dP{m;t) 
dt 



D 

w 

"o 



£[(' 



1)P(. 



l,n e + 1, 



;f) - riiP(. . . ,rii,n e 



;0 



(2.3) 
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where {e} denotes the set of nearest-neighbor sites adjacent to i, 



dP(rn;t) 



Of 



dP{m;t) 



(m + k){n i + k-l)...(n l + l)P(...,n l + k,...;t) 
rii(rii - 1) . . . (m - k + 1) P(. . . , m, . . . ; t) , 



dt 



h {k-\)d 

(m — m) P{. . . ,rii — m, . . . ; t) — m P(. . . ,rii, . . . ; t) 



and finally 



dP(m;t) 



dt 



/';....», • /';....//,....:/! 



(2.4) 
(2.5) 

(2.6) 



The precise form of the initial state is not important as we shall largely be considering stationary state properties. 
For convenience it is simplest to choose an imcorrelated Poisson distribution 



P({m};t) = e 



- p -Na(0) 



n 



n 



N A (0) 
N 



(2.7) 



As the configurations are given entirely in terms of the occupation numbers n^, this calls for a representation in 
terms of second-quantized bosonic operators 



whose effect on the state 



[ffli,ffl]-] = Sij , at |0) = , 



IW) = IlK t riQ) 



is to raise or lower the site occupation number by one, respectively: 

at | . . .m . . .)= I . . .rii + 1 . . .} , 
ai| . . .nj . ..)= rii | .. .rii — I. . .) . 

Then, upon defining the time-dependent state vector 



the master equation (2.2)— (2.6) may be recast into an 'imaginary-time' Schrodinger equation 



with the normal-ordered 'hamiltonian' 



D 



(2.8) 
(2.9) 



(2.10) 
(2.11) 



(2.12) 



(2.13) 



mwlh {«,}] = E H * = E -p E a J - a ^ 



Aa- 



,(k-l)d 



i-K T ) 



i - (4r 



ai en + t 



t\2 



1 - ( a i ) 



Using the initial state (2.7), the formal solution to Eq. (2.13) reads 

|$(t)) = er Ht |$(0)> = e~ NA ^ e~ Ht e"° E> a ' |0) . 



(2.14) 



(2.15) 



Our aim is of course to compute time-dependent expectation values of observables A, which may be defined in 
terms of the configuration probability according to 
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(A(t)) = Y l M{n i })P({n i };t) ; 

{m} 



for this, we need the projection state 



(PI 



n 



<P|0) = 



(2.16) 



(2.17) 



For then, as a consequence of (0| e ai being a left eigenstate of a\ with eigenvalue 1, one has (P\ Q({a\}, {a^}) = 
(P[Q({1}, {a>i}) for any normal-ordered polynomial Q of the ladder operators. Therefore, the expectation value 
( 2. 16[) may be written as 



(A(t)) = (P\ A({ ai }) |*(t)> = e-^W (0| A({oi}) eS, Q * e -*« e "° £, a ! | ) 
Notice also that probability conservation requires 

= <0| a * H({a\}, {aj) = (0| e^, * #({!}, {a J) , 



(2.18) 



(2.19) 



and therefore the hamiltonian generally has to vanish when the creation operators a\ are formally set to 1. We remark 
that in order to calculate 'inclusive' probabilities, such as the expectation value of the local density rij(t) = (a\(t)ai(t)), 
it is convenient to commute the factor e^-^ *^ through the hamiltonian in Eq. (2. IS), which is equivalent to a shift 



of all a\ — > 1 + aj, because then one just has to compute a vacuum expectation value of a normal-ordered product, 
for which Wick's theorem applies. Yet, this is not necessary for the evaluation of 'exclusive' quantities, as, e.g., the 
probability that site i is occupied by one particle, while all the other sites are empty, {S rii _i Y[j=a ^nj-,o}- Certainly this 
shift is not required for the renormalization of the model, and may in fact be even dange rous (as we shall see below). 

In order to arrive at a field-theoretic repre senta tion, we now view the computation of ( |2 . 1 8|) as a bosonic quantum 
many-particle problem with the hamiltonian (2.14), and again follow standard procedures pqjto write the expectation 
value as a coherent-state path integral 



(A(t )) 



/n i #,#,i({^})e- s '*^- ;t °' 



with the effective action 

S$i,ipi;to] = ^2 



dt 



&(*) ^Mt) + Hi{{lpi(t)}, {l/3i(t)}) 



- n o i>i(0) 



(2.20) 



(2.21) 



Finally, we perform the formal continuum limit, using the (assumed) inversion symmetry of the lattice, 

Y,^K d ld d x 7 ifc(t)-6gV(x,t), &(*)-► ^(x,t), 

i 

n 



b d n , ^[i(t)-^WH^ +2 V^(x,t) 



(2.22) 



which results in the action 



d d x 



dt \ ip(x,t) 



TT7 - DV 2 

at 



1-V(x,t) r 



^(x,i)-A fc [l-^(x,t) fc ]v(x,i) fe + 
$(x, t)i>{x, t) + r\l- $(x, t) 2 ^ 



-i/j(x,t ) - n o -0(x,O) 



(2.23) 



According to Eq. (2.20), S has to be dimcnsionless, and using Eq. (1.13) we find the scaling dimensions of the fields 
to be 



6 



$(x,t)] = «°, [^(x,t)] 



(2.24) 



and thus recover Eqs. (HTj) and ( 1.15 ). As mentioned in Sec. §. fluctuation effects are therefore expected to be 
important in physical dimensions in the cases k — 2 (d < 2) a nd k = 3 (d < 1). 

It is important to realize that in the above derivation of Eq. ( [2.23 ) , the form of the nonlinear terms followed directly 
from the master equation, and no further assumptions had to be made about these fluctuation contributions. On the 
other hand, the mean-field rate equation may be recovered from this action by considering the 'classical' equations 
of motion for the fields ip and ip, as given by the stationarity conditions 



ss 

Sip 

ss_ 

5ip 



d_ 

Bi 



DV Z ip ~ k\ k 



i - v> 



ip 



k-1 



l - V" 



d_ 

at 



DV' Z ) ip + k\ k ip k - 1 ip k + a 



l-(m + l)^ m V - 2tV3 - n 8(t) 



Obviously, Eq. (2.25) has the stationary homogeneous solution ip = 1, and inserting this into (2.26) yields 
£>V(x,f) 



dt 



= D V 2 V(x, t)-k\ k ip(x, tf +ma m V»(x, t) + 2r + n S(t) , 



(2.25) 
(2.26) 

(2.27) 



which, upon identifying ip(x,t) with a c oarse-grained local density (which is only possible on this mean-field level), 
is the natural generalization of Eq. (1.4). 

Notice that for even k and even m the hamiltonian part (in brackets) of the action (2.23), which does not explicitly 



depend on the final or initial state, is invariant under simultaneous 'parity' transformations of the fields, 

^(x, t) - -^(x, t) , V(x, t) -> -^(x, t) . (2.28) 

Physically, this corresponds to the fact that for even k and m the particle number is locally conserved modulo 2 
by both the annihilation and branching processes. However, if one expands the action (2.23) about the stationary 
solution according to 

^(x,t) = l + ^(x,t) , (2.29) 
corresponding to commuting the factor e^-'» a> through the hamiltonian in Eq. ( 2. IS ), one arrives at the new action 







jd d x 


JJdt(^{x,t) 







d_ 

at 



^ - DV 2 - ma„ 



^(x,t) + A fc ^Q 



V3(x,^V(x,i) fc 



m+l 



m + 1 
I 



—o- m 

i=i 

-V(x,0) -n o ^(x,0) 



■0(x, t) 1 t) - 2r -0(x, t)~r 4>(x, tf 



(2.30) 



in which the symmetry ( 2.28| ) for even k and m is obscured. Even worse, if in accordance with usual naive power 
counting arguments only the three-point vertices are kept and all the higher nonlinearities are disregarded, this hidden 
symmetry is lost completely. (Reasoning along these lines led a previous investigation to the erroneous conclusion 
that a field-theoretic treatment would predict that the dynamic phase transition for BARW with k = 2 should be 
in the directed percolation universality class irrespective of the parity of m [juj.) For a consistent renormalization 
group analysis, it is of course imperative to preserve all the symmetries of the problem. In the present case, this may 
be done by observing that the RG equations themselves (as opposed to the calculations of observables such as the 
density) should be independent of which basis is used, and it is therefore possible, and, indeed, necessary, to perform 
the computations in th e representation of the model in which the symmetry for even k and m is manifest, namely the 
unshifted action ( 2.23| ). 



III. BARW WITH TWO-PARTICLE ANNIHILATION NEAR TWO DIMENSIONS 
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A. RG eigenvalues at the annihilation fixed point and generation of new processes 



We now begin the inve stigation of fluctuation effects for BARW with two-particle annihilation (fc = 2, A2 = A), 
based on the action ( [2.23 ). For < t < t a , this field theory is characterized by a propagator Q(t) exp [— (Dq 2 + a)t\ 



(in the time domain) or (— iuj + Dq 2 + cr m ) _1 (in the frequency domain), and is graphically represented by a directed 
line, see Fig. |l](a), and by the vertices corresponding to the annihilation [Fig. |l](b,c)], branching [Fig. 0(d)], and 
pair production [Fig. 0(e)] reactions, respectively. With these elements, a systematic perturbation expansion may be 
constructed applying the usual techniques of quantum field theory [ p9[ . Taking into account fluctuation effects, we 
then define the renormalization constants Z\, Z„ m , and Z r according to 

£ = Z x XC d K d - 2 /D , (3.1) 
s m = Z am a m kT 2 jD , (3.2) 
T = Z T Tn- {d+2) /D . (3.3) 

The dimensionless renormalized couplings £ and s m introduced here are to be determined from the vertex functions 

r^^-q/ 2 ' " w / 2 ' -q/ 2 ' ~ w / 2 ; q/ 2 ' w / 2 ' i/ 2 > "V 2 ) and r f w(-q/ 2 ' ~ w / 2 ; _< i/ 2 > - u / 2 > 3c */ 2 ' 3cj / 2 )' 

respectively, at some specified sym metric normalization point setting a momentum scale k, e.g., <? 2 /4 = k 2 , uj = 0, 
or q = 0, iuj = 2Dk 2 . In Eq. (ful), C d = T(2 - d/2)/2 d - 1 n d / 2 denotes a d-dependent geometric factor (Ci = 1/2, 
C2 = l/27r). The renormalized pair production rate T may be obtained from the vertex function in the time 
domain, by introducing a lower cutoff i, and taking the Laplace transform with respect to ti at iu> = 2Dk 2 . Similarly, 
the renormalized annihilation and branching rates £ and s m may also be computed from the vertex functions 
and , respectively, in the time domain with an upper cutoff t f . 

Near the upper critical dimension d c — 2, we infer the scaling behavior by studying first the ultraviolet (UV) 
divergen ces of the perturbation theory, to be absorbed in the renormalization constants Z\, Z am , and Z T defined 
in Eqs. (|3.l|) -( |3~3|) . From the solution of the RG (Callan-Symanzik) equations for the correlation functions via the 
method of characteristics, we derive running couplings given by the differential RG flow equations 



dl 

ds m (l) 



'" (/) --ft(Z), *(0)=*, (3.4) 



dl 
dT(l) 

dl 



= ~Cs m (0 s m (l) , s m (0) = s m , (3.5) 
-(t(I)T(I), T(0)=T, (3.6) 



which yield the change of t, s m , and T under scale transformations k — > k e (the asymptotic regime being approached 



as I — y 00). Here, the beta and zeta functions follow from the Z factors in Eqs. (J3 . l|) — (3.3) according to 



Pi(l) = K^-X, (3.7) 
C„W = (3-8) 

OK <J m 

^ T {l) = K^-hx - . (3.9) 

OK T 

A scale-invariant asymptotic regime is then described by a zero £* of the beta function, 

ft(O = 0, (3.10) 
and the corresponding power laws are given by the anomalous dimensions (RG eigenvalues) 

y* m = -C m (n, (3.11) 

Vt = "(Tit) (3-12) 

at such a fixed point of the RG transformation, which is infrared-stable if f3' e (£*) > 0. 

The obvious first question to be addressed is whether the particle creation reactions with rates a m and r remain 
relevant perturbations to the pure annihilation process, i.e., whether the critical point remains at a c — t c = 0, once 
fluctuation effects are taken into account for d < 2 dimensions; and if so, what the corresponding critical exponents are. 
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To one-loop order, the only UV-divergent Feynman diagrams at d c = 2 which contribute to the renormalization of the 
annihilation rate A, the branching rate a m , and the pair production rate r are depicted in Fig. |^(a-c). These diagrams 
are identical with the lowest-order contributions at the annihilation fixed point, where a m = r = 0, provided the 
mass term oc a m in the propagator is set to zero. Upon using the dimensional regularization and minimal subtraction 
scheme, one then finds the following renormalization constants in d = 2 — e dimensions, 



Z x = l- 

Z„„„ = 1 - 



A C d K~ 



D e ' 
m(m + 1) A G d K~ 
2 D T 



Z T = 1 

D 



A C dK - e 



(3.13) 
(3.14) 
(3.15) 



The factor m(m + 1)/2 for Z„ m simply originates in the number of possibilities to select the external legs to the right 
of the bubble in F ig. B ( b). 

From Eqs. (|3 . 1[) , (p^|), and (pTq), we infer the beta function 



(3 t {£)=£{d-2 + £) 



(3.16) 



In the pure annihilation model (a m = r = 0), this is actually an exact result, which may be obtained either by 
summing a geometric series of the bubble diagrams in Fig. ^(a), or by writing down the equivalent Bethe-Salpeter 
equation for the annihilation vertex |J , leading to 



A C dK d - 2 
D 2-d 



and hence ( 3.16 ). For d > 2 the stable fixed point is just the 'Gaussian' fixed point I* 
non-trivial annihilation fixed point 



t 



2-d. 



T* = 



(3.17) 

0, while for d < 2 the 
(3.18) 



becomes stable and governs the asymptotic behavior. E.g., as there is no renormalization of the propagator, and hence 
neither field nor diffusion constant renormalization in the pure annihilation model, the solution of the RG equation 
for the average density may be written as 



n{n,£,t) = K d e~ dl n 



(tn 2 



-21 



1(1) 



(3.19) 



Near the fixed point (3.18), £(l) — > £*, we th en match e = tn , and arrive at the already cited exact result (1.16) 



Notice that in terms of the fixed point (3.18) the renormalization constant can be written as 

Z x = 1 - £/£* , 



(3.20) 



and h ence £ = £* corresponds to an infinite bare annihilation rate A, c.f. Eq. (3.1) 
(3.14), and ( |3 .15 ) yield the zeta functions 



(t(1) 



m(m + 1) 



£ + 0{£ 2 ) 



-(d + 2)+£ + 0(£ 2 ) 



Moreover, Eqs. (fj), (gj), 

(3.21) 
(3.22) 



and according to Eqs. (3.11) and ( 3.12[ ) the RG eigenvalues for the branching and pair production processes at the 
pure annihilation fixed point (3.1? 



become for d < 2 

m(m + 1) 



0(e 2 



y T =4-2e + 0(e 2 



(3.23) 
(3.24) 



Thus, both these couplings remain indeed relevant near d c = 2. However, it is important to note that in addition 
to the original processes all the lower branching reactions with m — 2, m — 4, . . . offspring particles become generated 
via fluctuations involving combinations of branching and annihilation processes, see Fig. 0(d). In a consistent RG 
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treatment all these additional reactions must be incorporated in the effective action from the start. In the case of odd 
m, the most important of these new processes is the single-particle annihilation A — ► 0. Upon assigning to that the 
decay rate n, 



decay : A — > , rate /x , 

which has the scaling dimension 

M = 

[compare Eq. ( [1.151 )1, the complete effective action reads 



(3.25) 
(3.26) 



m odd : S[ij>, ip; to] 



d d x 



dt\ i/>(x,i) 



(m+l)/2 



tp(x,t)-X l-^(x,i) 2 ip(x,t) 



l-i>{x,t) fl>(x,t)+ ^2 o-21-i l-V^x^) 2 '- 1 i)(x,t)i>(x,t) 



l=l 



1 - V3(x, tf - V>(x, t ) - n 4>{x, 0) 



(3.27) 



The emergence of the spontaneous decay ( |3.25 ) alters the behavior of the system drastically, and in fact leads to a shift 
of the critical point to a c > for d < 2, with the ensuing dynamic phase transition being in the directed-percolation 
universality class, and n(t) approaching its asympto tic va lue exponentially in both the active and ina ctiv e (absorbing) 
phases. We shall defer the discussion of the action ( 3.27 ) and the derivation of these results to Sec. VII . 

For even m, on the other hand, just the lower even-offspring branching processes are generated, and the effective 
action becomes 



even : S*^, -0; t ] = d 



dtl ^>(x,*) 



m/2 



ip(x,t)-X l-ip(x,t) 2 ^(x,t) 2 + 



1=1 



-ip(x,t ) - n o ip{x,0) 



(3.28) 



However, Eq. (3.23) implies that the most relevant of these branching reactions is actually the one with smallest m, 



i.e., m = 2, with rate oi = o and RG eigenvalue 



y a = 2 - 3e + 0(e 2 ) 



(3.29) 



The branching process with m = 2 will therefore describe the entire universality class of BARW with even offspring 
number m, clearly distinct from the universality class of BARW with odd m as a consequence of the underlying 
'parity' symmetry (2.28), i.e., the conservation of particle number modulo 2. In the following chapter, we shall derive 



the critical exponents for the BARW with m = 2 near two dimensions, and furthermore discuss the logarithmic 
corrections induced by the marginality of the annihilation vertex in d c = 2. 



B. The case m = 2: Critical exponents near d = 2 and logarithmic corrections 



We have just argued that the generic universality class of the dynamic phase transition for BARW with even 
offspring will be characterized by the critical exp onent s for m = 2. To one-loop order [see Figs. H(a-c) 1 , i.e., to first 
order in e = 2 — d, we found the RG eigenvalues ( 3.2E ) and ( 3.24 ) at the annihilation fixed point ( pTq ). Hence both 
a and r are relevant n ear d c = 2, and the critical point remains at its mean-field location a c = t c = 0. Upon utilizing 
the mean-field result ([[13) , we write the general solution of the RG equation for the density as 



n(n, £, s, T, t) — K d e~ 



sJV) 
£{l) 



i(t k 2 e~ 21 l{l) , s(l)/£(l) , T(l)/£(l) 



(3.30) 
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the pure annihilation result ( |3.19 ) is then recovered as s — > and T — > 0, if we demand that h(x, y, 0) — ► y 1 as 
y — > a nd n(a^) = lim^ofe/ ^(^M/j 0)]. Therefore at the critical point the density decays as n(t) oc t~ d / 2 , i.e., using 
Eq. ( 



d/2 



(3.31) 



Fo r non zero branching and pair production rate, but in the vicinity of the annihilation fixed point £(/) — > 
Eq. ( 3.30 ) reduces to 



n(K,i,8,T,t) (xse (v "- d)l h(tK 2 e~ 2l t, se y "\ Te Vrl 



As lengths scale as e l , the corr elati on leng th diverges according to £ oc s as s 
therefore we identify [see Eqs. ([Tl(]), (|L12|) ] 

v = \jy a , = l/y T . 

Finally, the matching condition e 21 = tn 2 leads to 

n{K,£, Sl T,t) oc st (v '- d)/2 h(t, st v " /2 , Tt v ^ 2 ^ , 



(3.32) 

and £ oc s^ 1 /^ as T -> 0; 

(3.33) 



which, upon comparison with ( p. 33] ) and (1.11) immediately implies 

z = 2 . 

Furthermore, f or T = and t — ► oo the t dependence has to cancel, i.e., h(£*,x, 0) oc as x 

the definition (1.9) we conclude that 

f3 = d/y a =dv = zvot . 



(3.34) 

(3.35) 
oo. Hence with 

(3.36) 



Notice that the relations ( |3.3l| ), fl3.33| ), ( |3.35| ) and ( |3.36[ ) are ezaci provided that cr c = t c = 0. Only the RG eigenvalues 
y a and y T need to be computed perturbati vely, e.g., in an expansion near d c = 2. 

Yet, assuming that the one-loop result ( |3.29| ) may be extrapolated with reasonable accuracy to some finite value 
of e = 2 — d, one finds that the branching reaction actually appears to become irrelevant for e > 2/3, i.e., below a 
supposedly new critical dimension 



d' c = 4/3 



(3.37) 



Notice also that the exponents v and diverge as d' c is approached from above. This would suggest the possibility 
of an entire inactive phase to emerge for d < d' c and a < a Cl i.e., the dynamic phase transition would be shifted to 
some positive critical value of the branching rate. For a < <j c the effective branching rate would then scale to zero 
under renormalization, and the power laws of the pure annihilation theory would apply. In such a situation the above 
scaling relations and values of the critical exponents would have to be modified considerably. 

In order to pursue this issue further, we now determine the RG eigenvalue y a and y T at the annihilation fixed 
point to second order in e. For that purpose, we need the UV-divergent two-loop diagrams for the annihilation and 
branching vertices, as depicted in Fig. |^(a,b). Taking the external momenta or frequencies at a symmetry point and 
using cith er of them to fix the normalization scale, one may readily compute the Z factors for A and a, as defined in 
Eqs. (|3.l[) and (3.2); using dimensional regularization and minimal subtraction, one eventually finds 



Zx 



D T 



Z CT =l-3- 



A 2 C 2 K~ 2 ' 
D 2 e 2 ' 

+ 6— d 

D 2 e 2 



— In — r~77 
2 3D 2 



(3.38) 
(3.39) 



Furthermore, by introducing a lower cutoff U for the diagrams contributing to the renormalization of the pair pro- 
duction rate, see Fig. ||(c), one finds 



A C d K~ 
D e 



X 2 C 2 k- 2 " 
D 2 e 2 



(3.40) 
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As to be expected, Eqs. (3.7) and ( 3.38 ) result in the same beta function ( |3.16D a s to one-loop order, for the 
diagrams in Fig. |^(a) are precisely those of the pure an nihil ation t heory , for which (|3.16| ) holds to all orders in £, and 
t her efore the fixed point is again £* — e. With Eqs. (3.39) and ( 3.40| ), the zeta functions for the branching vertex 
(3.8) and for the pair production rate (3.9) become 



2 + 3^-3 In - £ 2 + 0{£ 3 ) , 



CtW = -(d + 2) + e + o(f) 



(3.41) 
(3.42) 



and thus 



y a = 2-3e + 3 In- e 2 + 0{e 3 ) , (3.43) 
y T = 4-2e + 0(e 3 ) , (3.44) 

with t he RG eigenvalue for the pair production rate at the annihilation fixed point being identical with the one-loop 



result (3.24 



Eq. (3.43) for y a would indicate that the branching process becomes irrelevant for d < d" with 



1-Jl 



d" = 2 



8 i n i 

3 3 



2 In- 



1.1 



(3.45) 



and thus imply that in the interesting physical dimension d = 1 a stable inactive phase exists. Of course, an expansion 
near the upper critical dimension d c — 2 can by no means conclusively justify such an assertion. But we shall see in the 
following Sec. IV that a direct calculation in one dimension yields that surprisingly the one-loop result ( 3.23] ) becomes 
exact (with e = 1), and indeed the branching processes with even m are irrelevant in d = 1 near the annihilation fixed 
point. 

Yet, clearly in two dimensions all the previous results apply, and due to the marginality of the annihilation vertex 
induce logarithmic corrections to the mean-field critical exponents. With Eqs. (3.4) and ( 3.16| ) we find the running 
annihilation rate 



£{1) 



I + £{1)1 ' 

and upon inserting this into Eqs. (3J5), ( |3.6| ) and ( 3.21 ), ( 3.22| ) we may solve for the flows of s and T as well, 

=2i 



2 : 



s(l) 
T(l) 



s e 



[l + £(l)l] 

Te 4 ' 
1 + £{1)1 ' 



3 ' 



(3.46) 

(3.47) 
(3.48) 



We now employ Eq. ( 3.30| ) to infer the scaling behavior in the asymptotic limit I — > oo, where £{l) oc I 1 , s(l) oc e 21 1 3 , 
and T(l) oc e Al i -1 , and thus 



n{n, £,s,T,t)^ k 2 sT 2 n (t k 2 e~ 21 1~ 1 



5 2l r 2 , Te 41 



(3.49) 



Consequently, upon identifying £ with e l again, and by matching the second and third argument of n to one, 
respectively, we find 



T = : £ oc s- 1/2 ln(l/s) , 
s = : £ oc T" 1 / 4 ; 

in addition, from the prefactor of n we infer, using 21 ks ln(l/s), 

T = : n s (x s/[ln(l/s)] 2 



(3.50) 
(3.51) 



(3.52) 



Finally, at the critical point s = T = 0, Eqs. (3.49) or (3.19) become 
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n[K, £, t) = k 2 e~ 21 h(t n 2 e~ 21 r 1 ^ 
Matching the argument of h to a constant then leads to 



T = 



n(t) oc r 1 lnt 



(3.53) 



(3.54) 



which is consistent with the solution of the mean-field rate equation n = —\(t)n 2 with a time-dependent annihilation 
rate A(i) oc 1/lnt, sec Eq. (3.46) ||. For large initial densities no 3> n s , this logarithmic correction to the mean- 
field power-law relaxation may be observable even for nonzero branching rate, in an intermediate time window 



(2 A no)" 1 < t < (2 a)- 1 . Similarly, the non-trivial logarithmic corrections corresponding to Eq. J3.47| ) should 



appear in the opposite limit n -C n s for times (2er) 1 <C t <C (2\n ) 1 , inducing a subexponential increase 



n(t) oc exp 



2at 

(Lot) 5 



(3.55) 



At very long times, however, any positive branching rate will eventually lead to a purely exponential approach to 
n s . This crossover to an asymptotically Gaussian theory may actually be described by simply keeping the mass term 
a in the propagator for the one-loop diagrams depicted in Fig. |](a-c), which then results in the following beta and 
zeta functions, 



(3.56) 
(3.57) 
(3.58) 



p t (l,s) =tl d-2 + 
(s(l,s) = -2 



(1 + s) 2 - d / 2 
3£ 



(1 + S ) 2 - d /2 



( T (e lS ) = -(d + 2) + 



(1 + S ) 2 - d /2 



(for a related discussion of the crossover from a critical point to a Gaussian theory, see Ref. fl3 
coupling emerging here is 



The effective 
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(1 + , S )2-rf/2 



(3.59) 



and its flowing counterpart g[l) will vanish as I — ► oo because of the divergence of s(l), as long as s remains a relevant 
coupling. In this situation, i.e., for d > d' r , the Gaussian fixed point g* = describing the active phase is approached, 
and the above flow functions (3.56)— (3.58) assume their mean-field values, see Sec. |y|. 

In summary, we have found that near two dimensions the BARW with even offspring, described by the universality 
class of m = 2, are characterized by a critical point at vanishing branching and pair production rate, with critical 
exponents a = d/2, z — 2, and only two non-trivial exponents v and v T , which we have computed to second and first 
order in e = 2 — d, respectively. We have also calculated the leading terms for the logarithmic corrections appearing 
direct ly in two dimensions, and identified the active phase with a Gaussian fixed point in terms of the effective coupling 
( 3.59 ). However, there remains the possibility that in one dimension the branching process may actually be irrelevant 
near the annihilation fixed point, opening up an entire inactive phase for < a < a c . 



IV. BARW WITH TWO-PARTICLE ANNIHILATION AND EVEN NUMBER OF OFFSPRING IN ONE 

DIMENSION 

A. Computation of j/ CTm directly in d = 1 



In Sec. Ill it was shown that the RG eigenvalue y„ m which characterizes the relevance of the branching rate at 
(7 = may be computed within the e = 2 — d expansion, and, although it is relevant in d — 2, y am starts to decrease 
as e increases, indicating that it may eventually change sign. In this section we compute its value exactly in d = 1, 
by various methods, and show that this is indeed the case. 

In order to do this, we need to compute some dimensionless physical quantity, at some particular length or time 
scale characterized by the normalization scale n, as a function of the bare parameters a m and A, and understand 
how these latter should change as k is varied, keeping the physical quantities fixed. Since we are only int erested in 
the 0(s m ) terms in the RG equations (where s m denotes the dimensionless counterpart of a m , see Sec. Ill), we need 
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to compute these physical quantities only to first order in a m , although the computation should be carried out, in 
principle, to all orders in the annihilation rate A, since for d = 1 A is not small at the annihilation fixed point. Rather 



than the renormalized vertex functions considered in Sec. Ill, we shall consider a physically more accessible quantity, 
namely the probability P m+ i{t) that, given there is just one particle at t = 0, there are exactly m+ 1 particles found 
at time t. For sufficiently small t, we have P m +i(t) ~ o m t. Thus, we may define a renormalized version of a m in 
terms of the derivative dP m+ i(t) /dt, evaluated at some later time t = 1/Dk 2 . This is clearly quite different from the 
definition used in Sec. Ill; however, for universal quantities such as ycr m , it should yield identical results, 
make this renormalized branching rate dimensionless, we then define 



In order to 



t=l/Dn 2 



(4.1) 



where the tilde emphasizes the difference from the definition used in Sec. III. 

Pm+i{t) is given by the sum of all diagrams in the unshifted theory given by Eq. ( 2.23| ), which, in time-ordered 
perturbation theory, have exactly one particle at t = and m + 1 particles when cut at time t. To first order in 
cr m , there is only one branching vertex. Even so, the set of diagrams which may contribute to this may be very 
complicated. An important set is illustrated in Fig. ^. It is clear that, once the initial branching has taken place, 
the remaining portion of any such diagram corresponds to a solution of the quantum (m + l)-body problem for 
non-relativistic bosons interacting with a short-range repulsive potential of strength A. For general values of d and A 
this is intractable. However, in one dimension, and in the limit A — > oo, these bosons should behave like free fermions 
and therefore the problem be solvable. 

Fortunately, it is just this limit which is appropriate to discuss the behavior at the annihilation fixed point £ — £*, 
for this corresponds to a divergent bare coupling A. For d < 2 there are no ultraviolet divergences in the continuum 
theory, so that, on dimensional grounds, in one dimension P m +i(t) has the form (a m / A 2 ) F(X 2 t/ D) + 0(a^ n ), where 
F is some scaling function. Hence [c.f. Eqs. (3.2), (3.8)] Z am = F'(X 2 /D 2 k 2 ) + 0(s m ), and therefore 



1 



A 2 F"(X 2 /D 2 k 2 ) 



D 2 n 2 F / {X 2 /D 2 k 2 ) 



In this expression A is supposed to be expressed in terms of the renormalized coupling by 

A £*£ 



2Dn 



£ ' 



(4.2) 



(4.3) 



c.f. Eqs. (3.1) and ( 3.20 ). From these equations we see that, as £ — > £*, X — ► oo, and that this is the same as the limit 
t —> oo. 

This limit is more easily understood in the model defined on a lattice. For then two particles annihilate with 
probability one when they land on the same site. In this case, the branching process does not make sense unless 
the offspring are placed on different sites. Thus, we consider a slightly different version of the model in which the 
dimensionful parameter A is replaced by the lattice spacing bo- Once the m particles are placed on neighboring sites, 
they execute independent random walks until two or more of them fall on the same site and they annihilate. Let 
S m +i(t) be the probability that all m + 1 particles, initially placed on neighboring sites, survive until time t. Then 
we have the simple relation that 



P 



m+l 



(*) 



'm+l 



(i - t')dt' + 0{a 2 m ) , 



(4.4) 



so that s m = a m S m+ i(t = 1/Dn 2 ) + 0{a 2 m )- 

The asymptotic behavior of S m +i is simple to compute. This problem has been studied in d = 1 in several contexts 
before |3l]]. Denoting the coordinates of the particles by (x\,X2, ■ ■ ■ ,x m +i), we observe that this vector undergoes 
an isotropic random walk in the region Xj+\ — Xj > of an (m + l)-dimensional space, with absorbing boundaries 
along the hyperplanes a^+i = Xj, beginning at the point X® = jb . S m+ i(t) is just the survival probability of this 
walk. The Green's function for this problem is A m +i G^ m+1 \xi — x\, . . . , x m — x m ;i), where A m +\ is the completely 
antisymmetrizing operator on the m + l coordinates Xj , which ensures that the boundary condition is satisfied, and 
Q( m + 1 ) i s th e usual (m + l)-dimensional lattice Green's function on Z m+1 . In the long-time limit we may use a 
continuum approximation to this, with care: The survival probability is then given by integrating the full Green's 
function over the region xj + i > xj. The result must be a function of only the scaled variables —Xj)f(Dt) 1 / 2 , each 

of which is Oibo/iDt) 1 ^ 2 ) <C 1. The antisymmetry requires that the lowest term in an expansion in these variables 
must have the form 



14 



l<j 



This gives the result for finite lattice spacing, in the limit A — > oo. The corresponding result in the continuum limit for 
finite A is obtained by the simple replacement 60 — * (2/A). The precise coefficient (which is not, in any case, relevant 
to the present calculation) may be obtained by studying the case m = 2, which is solvable in both limits. 
Putting these results together, we then find that, in d = 1, 

y am = 2 - m(m + l)/2 , (4.6) 

exactly. 

The form of the result in Eq. ( |4.5| ) reflects the fermionic nature of the particles for infinite A. In fact, Eq. ( |4.6D may 
be derived more simply, if more formally, using this. Once the branching event has occurred, the particles propagate 
as an (m + l)-body fermionic state. If we express this in terms of anticommuting annihilation and creation operators 
C; and cj on each lattice site, the branching term in the hamiltonian has the form 



H b = a m y\ TT 4. . \ci (4.7) 




[for m even; in the case of to odd, the product runs from — (to + l)/2 to (m + l)/2, omitting j = 0.] The continuum 
limit of this expression is found by expanding in powers of jb$. This will be different from the bosonic case because 
the anticommuting nature of the c 1 allows each derivative to appear only once. Thus the lowest order, most relevant, 
term has the form 

b l+2+...+ m c f ( 9( ,t) ( 9 2 c f) _ _ _ ( 9 m c f) c f ( 4 _ 8 ) 

so that this term is multiplied by an effective expansio n pa rameter 6g m <r m , and its dimension is modified in a 



way corresponding precisely to the second term in Eq. (4.6) 



B. Results from an exactly solvable case 



Finally, we show how this result may, in the case m = 2, be derived from an exactly solvable lattice version of the 
problem. This is the model discussed by Takayasu and Tretyakov p7[ , in which a site is chosen at random, and, if 
there is a particle there, it is either moved to the left or right (with probability p), or two offspring are placed on 
the neighboring sites (with probability 1 — p). If the particles are interpreted as Ising domain walls, this model is 
a discrete time version of a kinetic Ising model solved by Droz, Racz, and Schmidt pM, with Glauber dynamics at 
effectively zero temperature, and Kawasaki dynamics at effectively infinite temperature. In fact this model is always 
in an inactive phase (except for p = 0), and does not exhibit the non-trivial transition to an active phase which is one 
of our main concerns in this paper. As shown in Ref . p2| , it is necessary to include another parameter in the model, 
which enhances the branching rate relative to the diffusion process, in order to find such a transition. Nevertheless, 
for the purposes of extracting the universal eigenvalue y a , at the pure annihilation fixed point, this deficiency is 
unimportant. 

Takayasu and Tretyakov [jl7|] consider the probability Q r (t) that a given interval of length r contains an even number 
of particles. Since they assume translationally invariant initial conditions, it does not matter which interval. In the 
Ising spin language, this is simply related to the correlation function (s, s J+r ). They show that this satisfies the linear 
system of equations 



Q r (t + 1) - Q r (t) = (1/N) p - p)(Q r+1 - 2Q r + Q r+1 ) 
Qi(t + 1) - Qx{t) = (1/N) \(2 - p)Q 2 - 2Qi 



(r>2) 



(4.9) 
(4.10) 



where N is the number of sites. Although Takayasu and Tretyakov did not analyze this equation in detail, it is not 
difficult to take a continuum limit and extract the scaling behavior of the parameter 1 — p, which is the analog of 
the branching rate a = 02- Introduce Qo so that the second equation above has the same form as the first, extended 
down to r = 1. This means taking 



2(l-p)(Qi-Q )=pQ . 



(4.11) 
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On rescaling t — > t = t/N and r — > f — b$r, the continuum limit 1/N — > 0, f>o — * now gives a simple diffusion 
equation 



for all real r > 0, together with boundary condition 

2(1 -p)b ^=pQ (4.13) 

at r = 0. This corresponds then to a heat diffusion problem on the half-line r > 0, with a radiating boundary at 
r = 0. There is a length scale Lf, = (1 — p)bo implicit in the boundary condition. For r ^> the solution at late 
times will approach that with the Dirichlet condition Q(r = 0) = 0. However this will be modified for r <C Lt- The 



fact that 1 — p ~ a scales with length in this way immediately implies that y a — —1, consistent with Eq. (4.6). In 
fact, for this continuum model, this is true to all orders in a. This is because the non-trivial fixed point, determining 
the transition at p — 0, is infinitely far away in this scheme. 

Within this model, it is possible to determine P$(t) and hence check the calculation given earlier in this section. 
If there is initially only one particle in the system, then Q r (t — 0) = r/N . If p = 1 (no branching), this is an exact 
steady state solution, corresponding to a uniform heat current. Let us find the solution in perturbation theory in 
1 — p ~ a, by writing Q r (t) = (r/N) + (1 — p)f(r, t) + 0((1 — p) 2 )- Then / obeys the diffusion equation, vanishes at 
t = 0, and satisfies 2(1 — p)/N + 0((l — p) 2 ) = (1 — p)f at r = 0. It is thus given in terms of the Green's function 
Gn(t — t';r,r'), satisfying Dirichlet boundary conditions at r = 0, by 

/M) .j.fw-WV-jgfio-Lp,.-*^* . (4.14, 

To first order in u, there are only one or three particles in the system, so that the density is Qi{t) = (-Pi + 3Ps)/N — 
(1 + 2P 3 )/N. Hence P 3 (t) = (N/2)(l-p)f(r = b ,t) + C((l -p) 2 ), and the 0(1 - p) contribution to P^t) therefore 
scales as i~ 3 / 2 . Following through the RG procedure described earlier then leads to y a2 = — 1 as before. 

Within this model, it is possible to investigate the nature of the long-time behavior to all orders in ui ~ 1 — p. If 
p < 1, the boundary temperature Qo(t) will increase to a point where the heat current crossing the boundary equals 
that coming from large r. Thus Qo(t) ~ 2(1 — p)/N, so that the density of particles, Qx(t) ~ [1 + 2(1 — p)/p]/N. 
There is therefore a finite average number of particles in the whole system, which diverges as p — > 0. It is easy to 
check that this is the asymptotic solution for any finite odd initial number of particles. On the other hand, if this 
number is even, then initially Q r (0) — r(N — r)/N . In this case, the only steady state solution is Q r = 0, and it is 
straightforward to check that the particle density Qi(t) ~ t^ 1 / 2 , as expected in the inactive phase. 



V. TWO-OFFSPRING BARW: TRUNCATED LOOP EXPANSION 



A. RG flows to one— loop order 



In Sec. [II B we saw that near d c = 2 dimensions, the branching rate a for BARW with two offspring particles 
remains a relevant perturbation. Thus the critical point is at a c = 0, as in mean-field theory, and for any positive 
branching rate there exists only the active phase, which is governed by a Gaussian fixed point. On the other hand, 
the direct analysis of the one-dimensional model in Sec. E\l unambiguously established the existence of an inactive 
phase described by the power laws of the pure annihilation model. Hence, in the inactive phase, a must constitute 
an irrelevant operator in the RG sense. Clearly, this behavior cannot be constructed within an expansion about the 
upper critical dimension d c — 2. Yet, by employing a truncated loop expansion at fixed dimension we can establish 
a unifying RG framework for the above results. To one-loop order, we shall indeed find a new critical dimension 
d' c = 4/3, below which the inactive phase emerges. 

The complete set of tree and one-loop diagrams for the two- and four-point vertex functions is displayed in Fig. |^. 



As opposed to the analysis in Sec. [II, we now retain the full dependence of both the vertices and the propagators on 
the branching rate a. Upon applying an upper cutoff tf in time for the processes depicted in Figs. ^|(a,b), taking the 
Lapl ace tr ansf orm with respect to tf at the normalization point iuj = 2Dk 2 , and employing our previous definitions 



(3.1) and (3.2), one readily arrives at the renormalization constants 
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Zx = l- 
Z„ = 1 - 3 



X/D 



2-d (Ki + tT/D) 1 -*' 2 

_Cd A/D 

2-d 



{n 2 + a/D) 



l-d/2 



Consequently, with Eqs. (3.7), (3.8) and & = we find the zeta functions 



Ci{t,s) = d-2 
( s (£,s) = -2 + 



(l + s )a-d/2 
31 

(1 + s)2-d/2 ' 



(5.1) 
(5.2) 

(5.3) 
(5.4) 



For s = 0, of course, these reduce to Eqs. ( 3.16Q and (3.21) (for m = 2) as computed at the annihilation fixed point 
(3.18). Within this truncated one-loop expansion, we find for the RG eigenvalue of a at this fixed point 



Va = 3d - 4 



(5.5) 



which indicates that indeed for d < d' c = 4/3 the branching processes become irrelevant. Note that in one dimension, 
this happens to coincide with the exact result (4.6) for m = 2. For d > d' c , however, the annihilation fixed point is 
unstable with respect to branching, and s(l) will diverge as I — » oo. In this case we expect the flow of I and s to be 
described by their naive scaling dimensions. Indee d, w ith i l l) ~ e -{d-2)i an j gffl ^ e 2i ^ we nave £(},)/ s{l) 2 ~ d / 2 ~ 
e~ 2 ' — » 0, and the fluctuation contributions to Eqs. (5.3), (5.4) vanish. In Sec. Ill B , we had already anticipated these 
flow functions describing the crossover to the asymptotically Gaussian model that characterizes the active phase. 

Before discussing the ensuing flow equations in more detail, we compute the renormalization of the pair production 
rate as well. Notice that in this calculation at fixed dimension, we cannot just restrict ourselves to those diagrams 
that become ultraviolet divergent at d c = 2. The three graphs in Fig. ||(c), computed with a lower cutoff U, lead to 
the Z factor 



Z T = 1- 



c d 

2-d 



X/D 



3C d <jX 



(k 2 + a/D) 



!- d / 2 2-d D 2 K 2 



1 



1 



(k 2 + <r/D) 1 - d / 2 (a/D) 1 ^/ 2 



(5.6) 



For cr = 0, this obviously reduces to Eq. ( 3.15 ); on the other hand, as a 
brackets with the result 



oo, we may expand the term in square 



<a/D 



1 



C d 
2-d 



X/D 



X/D 



(k 2 + a/D) 1 - d/2 



3C d 

2 {a/Df- d l 2 



O 



X/D 



{a ID) 



2-d/2 



With the definition (3.9) this leads to 



C T (£,s) = -(d+2) + 



(l + s )2-d/2 ' 



(5.7) 



(5.8) 



which is identical to Eq. ( 3.58 ). Note that Eq. (5.8), which has been obtained to leading order in an expansion for 
large <j/Dk 2 , correctly incorporates both the inactive (s = 0) and active (s — ► oo) phases, as do Eqs. (5.2) and (|5.-j|). 

It is instructive to realize that one arrives at the same results for 
Q and ( s by considering the four-point vertex functions r^^^^(— q/2, — uj/2; — q/2, — ui/2; q/2, u/2; q/2, uj/2) and 
r^^(— q/2, —uj/2; —q/2, —uj/2; —q/2, — ui/2; 3q/2, 3w/2), respectively, at either the normalization point g 2 /4 = k 2 , 
lo = 0, or at q = 0, iui = 2Dk 2 . This comes about beca use the expressions from the triangular loops in Figs. ||(d,e) 
may be written in the form of the last term of Eq. (5.6), and hence do neither contribute at the annihilation nor at 
the Gaussian active fixed point. Finally, the loop in the fourth diagram in Fig. |^(d) carries no external momentum 
or frequency, and thus has no influence on the zeta function. As in the one-loop e expansion near d c = 2, there is 
neither field renormalization nor renormalization of the diffusion constant at this level, because the loop diagram for 
r,f^(q, u>) (with no upper cutoff in the time domain), see Fig. ||(b), carries no wave vector or frequency dependence. 

Returning to Eqs. (5.3), (5.4), and ( |5.8| ), it is clear that the relevant effective coupling is indeed g = £/(l + s) 2 ~ d / 2 , 
as defined in Eq. (3.5S). In the prospective inactive phase, with s — > 0, we have g —> £, the coupling of the pure 



annihilation model. On the other hand, if s diverges, we may construct the beta function for g as follows: 



PM^9 &-(2-d/2)& 



10 - id 



(5.9) 
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Thus, as explained above, despite both 1(1) and s(l) diverging according to their naive scaling dimensions, there is a 
stable Gaussian fixed point g* = describing the active phase. Yet, Eq. (5.S) yields the additional non-trivial fixed 
point 



■ 9 * = io^3d ■ (5 - 10) 

As 

dRJn\ 

= 2- (10 -3d) g (5.11) 



dg 



assumes the value —2 at g — g* , this new fixed point turns out to be infrared-unstable. As the fixed-point values for 
the new effective coupling arc limited by those of the annihilation rate, which in turn has as upper bound £* = 2 — d 



(for according to Eq. (4.3) this already corresponds to an infinitely large bare coupling), g* enters the physical regime 
for d < d' c — 4/3, precisely when the annihilation fixed point and the associated inactive phase become stable. The 
non-trivial unstable fixed point then governs the dynamical phase transition between the power-law inactive and the 
Gaussian active phases. 

The flow diagram in d = 1, as obtained by numerically solving the flow equations, is shown in Fig. ||. The flows 
in the upper left half describe the active phase with both s(l) — > oo and t(l) — > oo, but the effective coupling g 
approaching a Gaussian fixed point, g(l) — > 0; on the other hand, the trajectories in the lower right part of the figure 
correspond to the inactive phase, with s(l) — > and g(l) — > 1(1) — > £* — 2 — d (annihilation fixed point) as I — > oo. 
Notice the remarkable flows in the inactive phase, with the RG trajectories curling about the annihilation fixed point. 
The unstable fixed point with £(l) oc s(l) 2 ~ d / 2 (for large £ and s) appears as a separatrix separating the basins of 
attraction of the fixed points associated with the inactive and active phase, respectively. For any dimension d < d' cl 
the flow diagram would look qualitatively similar. At d = d' c = 4/3 (to this order) itself, there appears an attractive 
line of fixed points with Ce — Cs — 0, starting at the point £ = 2/3, s — 0. For d > d' c , all the flow trajectories tend to 
the Gaussian fixed point g* — (with both t, s — > oo). 

The ensuing phase diagram as function of space dimension is displayed in Fig. |?]. For d > d' c and any nonzero 
branching rate there exists only the active phase. For d > d c = 2, the transition at a c = is described by the mean- 
field exponents, see Sec. \ , whi le for d' c < d < 2 one finds the non-trivial critical exponents (logarithmic corrections at 



d c = 2) evaluated in Sec. [II B. Only for d < d' c , when g* < £* = 2 — d and the annihilation fixed point becomes stable 
for sufficiently small values of s, the inactive phase which is governed by the power laws of the pure annihilation model 
emerges. The critical behavior at the dynamic phase transition with a c > is characterized by exponents belonging 
to a novel universality class. In the following subsection, we shall analyze the scaling behavior at this transition, and 
compute the critical exponents to one-loop order. Although the actual values we find for these exponents turn out 
to be rather poor estimates, as a consequence of the absence of any small expansion parameter, we believe that the 
qualitative features of the phase diagram, and specifically the mechanism for how the inactive phase becomes possible, 
are correctly encoded in the truncated loop expansion. 



B. Scaling analysis and critical exponents 

We now w ant t o explore the critical behavior in the vicinity of the transition described by the unstable fixed point 
g = g*, Eq. ( 5.10| ). We thus introduce 



.9' - 9 



(5.12) 



and note that because of g*de(l)/dl — (3 g (g(l)) = —g*e(l) dj3 g (g) / dg\ _ „ , we identify the corresponding zeta function 



g=g 



dP g (g) 



dg 



(5.13) 



Collecting the results (|5.3|), (5.4), ( |5.8| ), and (5.11), we then find the following RG eigenvalues at the unstable fixed 
point g*, 
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He 

yx 



-C 



-Q = 

Va = -C = 
Vr = -Ct = 



(4-d)(4-3d) 

10 -3d 
2(4 - 3d) 
10 -3d ' 
16 + 4d - 3d 2 
10 -3d 



(5.14) 
(5.15) 

(5.16) 

(5.17) 



We remark that y\ and y a are not independent, but fixed by the condition that f3 g (g* > 0) = 0, and hence y\ — 
(2 — d/2)y a . Notice also that to one-loop order there are no anomalous dimensions for the diffusion constant or the 
fields themselves. 

In identifying these RG eigenvalues with the critical exponents, however, we have to be careful as now the critical 
point is shifted away from a c — to g c = g*, and we therefore cannot apply Eqs. ( |3.33 ) and (3.36). In general, the 
critical exponents will rather depend on y £ , as well as on y\ and y a . Only the exponent v T is not affected by these 
modifications, and therefore 



v r = l/y T 



(5.18) 



remains valid. In order to compute the exponents at r = 0, we write the solution of the RG equation for the density 
in the form 



n(n, £, s, t) = n d e~ dl h (k, £(l), s(l),t(n 2 /D)e- 21 ) 



(5.19) 



By itself, this is not sufficient to determine the exponents, since s(l) and £{l) do not themselves flow to fixed points 
and the dependence of n on these, even at early times, is not known a priori. In the active phase, however, we k now 
that they do flow to a region in which mean-field theory is valid. This predicts that the right hand side of (5.19) has 
the form s(l)/£{l) times a function of the combination t{n 2 / 1 D)e~ 21 s(l) . In the vicinity of g*, then, (5.19) becomes 



n(n, £, s, t) = n d e-^+y^y^ 1 jn(ee y ' 1 ,t K 2 e^ 2 " 2 ^ 1 s) 



(5.20) 



Upon choosing the matching condition e ' = lel 1 ^ 6 , this implies that we may identify the critical exponents, as 
defined in Sec. |, as 



= (d + yx - y<r)/y s 

v=(2-y a )/2y e , 
z = 2 . 



(5.21) 
(5.22) 
(5.23) 



Notice that lengths now scale as e {1 - y ° /2}l . The mean-field value for the dynamic exponent z naturally follows from 
the absence of diffusion constant renormalization. 

We emphasize that as a consequence of the appearance of the dangerous irrelevant variable 1/s, we cannot unam- 
biguously extract the exponents defining the power laws at the critical point itself, such as a. In addition, this also 
precludes us to provide a sound foundation based on the renormalization group for scaling relations like (3 — zvet 
]Tl1 . [r9|j22| for the dynamic phase transition for d < d' c . Such scaling relations would follow only if the mean-field 
dependence s(l)/£(l) were also to be valid as both s(l) and £(l) grow without bound along the critical line g = g* . 
If the fixed point occurred at finite values of these quantities, it would be reasonable to assume that the right hand 
side of Eq. ( |5.19| ) does have this mean-field form, as it is to be evaluated at some finite rescaled time. It would be 
interesting to look care fully fo r pos sible violations of these scaling laws within simulations. 

Upon inserting Eqs. ( 5.14 )-( p.l7 ), we find the following results as function of dimension, 



10 -3d 



10 -3d 



10- 3d 



16 + 4d - 3d 2 



, z = 2 



(5.24) 



At the borderline dimension d' c — 4/3, the v alue for the critical exponent v T = 3/8 coincides with the one found in 
the e expansion about d c = 2 (see Sec. |III B ). On the other hand, upon approaching d' c from below, one has v = 1/2 
and (3 = 2/3, which has to be contrasted with v, /3 — * oo as d [ d' c . In the more interesting physical dimension d = 1 
we find to one-loop order 



d=l : ^=3/7,^ = 2/7,^ = 17/7, 



(5.25) 
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and thus 



d=l : = 4/7, v = 3/7, v r = 7/17, z = 2 . (5.26) 



These numerical values, compared to the actual simulation results ]1 1| , |17| , |22| , |26[ are generally rather poor, with the 
remarkable exception of v r . Keeping in mind that the truncated loop expansion, with no small parameter at hand, 
constitutes an uncontrolled approximation scheme, the unsatisfying accuracy may not be too surprising. 

C. Some remarks regarding higher orders in the truncated loop expansion 

It is by no means clear that including higher loop orders in the above analysis would yield considerably better 
results for the critical exponents, given that there is no a priori small expansion parameter present at all. Yet some 
general remarks about a possible extension, say, to two-loop order, are in place here. 

To two-loop order, one would obviously expect both field renormalization and diffusion constant renormalization 
to appear, e.g., from the frequency and wave vector dependence of the loop in the fourth diagram of Fig. ^(b), or 
more generally, from the set of graphs depicted in Fig. | . Defining 

C, = «^mf, (5.27) 
Co^Afc^, (5.28) 



this will imply two additional RG eigenvalues y^p = — Q, and yo = —Cd- As m ay be inferred from Eqs. (5.19) and 
fl5.20 ), their appearance modifies the previous scaling relations ( |5.21 )-( 5.23| ) to 



(i = {d + y^ + yx-Va)/Ve , (5.29) 
u=(2 + y D -y a )/2y e , (5.30) 

-^4. (5-31) 
2 + yD - y<j 

and would thus allows for a non-trivial value of the dynamic exponent as well. 

It is, however, far from obvious that an extension of the above one-loop analysis to higher orders in the perturbation 
expansion is feasible. Certainly, at the annihilation fixed point describing the inactive phase, a loop expansion is well- 
defined to any order. Also, the Gaussian fixed point is trivially described by the naive scaling dimensions of the model 
parameters. Yet the unstable fixed point governing the dynamic phase transition requires a subtle balancing of the 
divergences of the couplings s(l) and 1(1), as encoded in the appearance of the new effective coupling g(l). To one-loop 
order, the expansion parameter in the bare theory is (\/D)/(k 2 + a / 'D) 1 ~ d / 2 , and hence the effective coupling in the 
flow function becomes g — 1/(1 + s) 2 ~ d / 2 . The two-loop diagrams, then, are proportional to (X/D) 2 /(k 2 + o/D) 2 ~ d , 
and thus their contributions to the zeta functions proportional to £ 2 /(l + s) 3 ~ d — g 2 (l + s) — which diverges at 
the transition. It appears, therefore, that in an expansion in 1/s, the leading terms would have to cancel in order 
to render the loop expansion well-defined. However, although it may be shown for Q that this is fact true for those 
two-loop diagrams that in the frequency domain are product of the one-loop terms, this appears not to be the case 
for the nested loop diagrams such as the fourth graph in Fig. |^(b), and perhaps the summation of an entire series of 
such diagrams, see Fig. |], might be indispensable. Despite some effort, we therefore did not succeed in extending our 
one-loop analysis further, or even to demonstrate that there exists a meaningful truncated two-loop theory at all. 

VI. GENERALIZATION TO N SPECIES OF PARTICLES 

We now generalize the BARW with k = m = 2 to N species of particles a = 1, . . . , N, with only identical particles 
annihilating each other, and two different types of branching processes, 

annihilation : A a + A a — > , rate A , (6.1) 

branching : A a — » A a + A a + A a , rate a , (6.2) 

A a — * A a + A 13 + A 13 (P^a), rate a'/(N - 1) , (6.3) 

pair production : — > A a + A a , rate r . (6-4) 
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The motivation behind such a generalization is that we found in Sec. |y] that in the interesting physical case of N = 1 
in d = 1 the loop expansion was uncontrolled and does not lead to reliable values for the exponents, even though 
it correctly predicts the existence of a transition. As in the case of iV-component magnets, then, one might hope 
that there exists an TV-species generalization of the model with m = 2 introduced in Sec. EL which becomes exactly 
solvable in the limit N — ► oo, yet which still retains the essential phys ics. Unfortunately this does not seem to be 
so. There appears to be no simple way of generalizing the action ( 3.28 ) (with m — 2) so that it possesses an O(N) 
symmetry as in the magnetic case. The highest symmetry that can be embedded while at the same time retaining the 
conservation modulo 2 of each species is Sn, the permutation symmetry of the N species. This is clearly exhibited 
by the allowe d re action s ab ove. Note, however, that the permutation symmetry allows different possible rates for the 
reactions in (6.2) and (6.3). The symmetries of this model are most clearly seen in the form of hamiltonian density 
in the second-quantized formalism: 



4>l 



1 VaV'c* 



N-l ^ 



\4>a - 1) i>p^p 



(6.5) 



When a' = (N — l)er, the branching terms do in fact exhibit an O(N) symmetry which generalizes the Z2 symmetry 
of the original model. However, the annihilation terms inevitably violate it, and any attempt to rectify this falls foul 
of the requirement of probability conservation. This has the consequence that the symmetry of the branching terms 
is not preserved under renormalization. Physically, this is because after the first process (6.2), any pair of the three 
outgoing particles may ann ihil ate, so the chances of observing three particles after some amount of time are reduced 
as compared with process ( |6.3[ ), where only one pair may annihilate. 

The two-loop calculation of Sec. IIIB is readily generalized to ./V particle 'flavors' by modifying the appropriate 
combinatorial facto rs in th e evaluation of the Feynman diagrams depicted in Fig. ||. Near two dime nsio ns, this yields 



the Z factors ( 3.38 ), ( 3.39 ) as before, while the graphs in Fig. ||(b) for the new branching process (O) give 



1 - — 



A_ C d K~ 
D e 



*1 
D 2 



(6.6) 



Consequently, with the definitions used in Sec. Ill, the beta function (3.16) and fixed point ( 3.18| ), as well as Eq. (3.41) 
remain unaltered, while the additional zeta function reads 



and therefore the RG eigenvalue of the branching rate a' at the annihilation fixed point is 

y a , =2-e + 0(e 3 ) . 



(6.7) 
(6.8) 

Upon comparing this result with ( |3.43] ), we see that the branching with identical offspring ( |6.2| ) is irrelevant as 
compared to the production of different particle 'flavors' (|6.3[): for I — ► 00, one has 



s(l)/s'(l) ~ exp{(y a - y a ,)l} = exp{[-2e + 31n(4/3)e 2 ] 1} 







(6.9) 



Therefore the process (6.3) is going to dominate the long-time behavior, and we may effectively set a = 0, but retain 



a 1 > 0, for any N > 1. The remaining Feynman diagrams are then identical to those that would be obtained in the 
limit N — > 00, and one therefore expects that actually all reactions wit h N > 1 are asymptotically described by the 
model with infinitely many particle species and sole branching process (6.3). 

The above model in the limit a = 0, or equivalently, for N — > 00, has a considerably simpler structure as compared to 
the N — 1 case, because 'nested' diagrams as, e.g., the last graph in Fig. ||(b) do not appear any more. Consequently, 
all fluctuation contributions to the propagator may be absorbed into a renormalization of the 'mass' a. Therefore one is 
left with the task of summing an infinite scries of branching 'bubble' diagrams ('cactus' graphs). The most convenient 
way of achieving this is to write down self-consistent Bethe-Salpeter type equations for the annihilation, branching, 
and pair production rates. Graphically, these are represented in Fig. |^(a-c), respectively. The corresponding analytic 
expressions read 



At? = --1^(0,0) = A 



^ = r^(o,o) = a' 
1 



1 - -E- 



1 



D J p P 2 + °r/D\ 

x r r 1 



TR 



r«(o,o) = A 



2 



D J p P 2 + o' R /D 
, X R f 1 



D J p P 2 +a' R /D 



(6.10) 
(6.11) 
(6.12) 
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which are readily solved by 



X R /X = <t'r/<t' = tr/t 



X 

1 + - 



1 



DJ P P 2 + v'r/D\ 

If we now evaluate the integral at the normalization point a' R /D = re 2 , we arrive at the exact result 

which immediately implies that Eqs. ( |3.16| ), ( 3.22 ), and (|6.7| ) hold to all orders in i, 

t (e) = £(d - 2 + i) , ( s ,(£) = -2 + £, C T (£)=-(d + 2)+e . 
Thus in this model, a' c = t c = 0, and the only non-trivial RG eigenvalues are 



y<r> 



Vr 



2d 



(6.13) 



(6.14) 



(6.15) 



(6.16) 



which determine the divergence of the correlation length as either a' —> or r — » 0. 

We have thereby demonstrated that the above iV-species generalization of the BARW with m = 2 leads to a new 
universality class, which for all N > 1 is asymptotically characterized by the model with N — > oo. In this model, the 
critical point remains at vanishing branching and pair production rate, and for d < 2 its critical behavior is governed 
by the exponents 



v=l/d, v T = l/2d, 2 = 2, a = d/2, (3 = 1 



(6.17) 



In d c = 2 dimensions, one finds a logarithmic correction for the de nsit y decay at the critical point, just as in Eq. (|3J54l) . 

The N — * oo limit of the model defined by the processes (6.1)-([0|) may in fact be solved directly, without 
going through the renormalization group machinery. This is based on the simple observation that in the reaction 
A a — > A a + A 13 + A 13 , the two products A 13 are far more likely to annihilate against each other, since they have the 
same flavor index, than against the remaining A a , or indeed any other particle in the system, which is unlikely to have 
a matching flavor index. In the limit N — > oo, this becomes overwhelmingly the case. Now the average particle density 
can change by branching processes, which will always be proportional to the existing local density, and annihilation 
processes. These latter are either between particles which are siblings, that is they are products of the same branching 
event in the past, and which will not, in this limit, annihilate against any other particle; or they are between particles 
which are from independent branching events. In this case they are unlikely to have the same flavor index, so that 
this term is suppressed by a factor 1/-ZV; however, it is proportional to the square of the mean density, which is itself 
O(N), and therefore it also gains a factor of N relative to the linear terms. If the particles come from independent 
branching events, we may also, to leading order in 1/N, neglect the correlation between these events. 

One may therefore write down an integral equation for the time evolution of the mean density n(t): 



dn 
~db 



2a n{t) - 2a [ L(t - t') n(t') dt' - X [ [1 - L(t - t')] n[t'f dt' , 
Jo Jo 



(6.18) 



where L(t — t')dt is the probability that two particles, created at time t' , annihilate in the interval (t,t + dt). Note 
that L{t) = —S' 2 {t), where, as in Sec. Saft — t') is the survival probability at time t for two isolated particles 
produced at time t' . The factor of 1 — L(t — t') in the last term ensures that the particles have not annihilated at 
some earlier time. 



Eq. (6.18) may be solved simply by Laplace transform. Defining h(p) = f °° n(t)e pt dt, and similarly for L(t), one 
finds 



pn(p) — n(0) = 1 — L(p) 2a'h{p) — Xn 2 (p) 

where L(p) is given by a simple sum of bubble diagrams as 

2XI( P ) 



Hp) 



l + 2XI{p) ' 



(6.19) 



(6.20) 



and I(p) = J. l/(p+ 2Dk 2 ). Although (6.18) cannot be solved in closed form, the various critical exponents may 
easily be deduced. If we begin from an initial state when n(0) is small, we may initially neglect the nonlinear term to 
find that, whenever a' > 0, h(p) has a pole for positive real p at the solution p = Po(a') of 
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p [1 + 2X1 (p)] = 2a' . (6.21) 

The density therefore increases exponentially at a rate e pot . When d > 2, I(p) is analytic at p = and we see 
that pa oc a' . When d < 2, po oc o' zv , where zv = 2/d, and for d — 2, po ~ ct'/ ln(l/cr'), all in accord with the 
renormalization group analysis presented above. For later times, this exponential growth is capped by the nonlinear 
term an d th e system reaches a steady state. Since the factor 1 — L(p) is common to both terms on the right hand 
side of (6.19), the steady state density behaves simply as cr /5 , where (3 = 1 irrespective of the value of d. 



VII. BARW WITH ODD NUMBER OF OFFSPRING 

A. Mapping to Reggeon field theory for directed percolation 

In this section, we consider the case of to odd, which turns out to be quite different from the case of even to. For 
odd m, there is no 'parity' symmetry and one might expect the transition from an inactive state, provided it occurs 
at some non-trivial value of the branching rate o~ m , to be in the directed percolation (DP) universality class, as a 
consequence of the generated particle decay processes A — > 0. This in fact is what we find in low dimensions, but 
there are subtle issues connected with the fact that these decay processes are themselves fluctuation-induced, and 
thus proportional to the branching rate o~ m ■ Therefore the very existence of a non-trivial transition must depend on 
fluctuation effects, which are expected to be important only when d < 2. As pointed out in Sec. |, the mean-field 
equation always predicts an active state for all a m > 0, and is expected to be valid, at least in the pure annihilation 
problem, for d > 2. Indeed we find that only for dimensions d < 2, fluctuations are strong enough to produce a nonzero 
critical branching rate a c > 0, and consequently a DP phase transition separating the active from the inactive phase. 
This is despite the fact that the upper critical dimension for the DP transition itself is d c = 4. 



For these reasons we investigate this problem in the vicinity of two dimensions. As was shown in Sec. [II, from 
an initial model with a given value of to > 1, under renormalization the other processes with numbers of offspring 
to — 2, to — 4, . ..,1,-1 are generated, compare Fig. ||(d). Therefore we shall begin by considering the si mplest case 



to = 1 . We start by writing the interaction terms appearing under the integral in the effective action Eq. ( 3.27| ) 



Hi„t = — A (i — ^ 2 ) ip 2 - <ti - 1) 4>ij - m(i - $)1> • (7.1) 

The last term corresponds to the process A — > 0, which, although not present in the original model, is generated by 
the combined processes A — > 2A, 2A — > 0. 

It is also important to record the form of the interaction in terms of the shifted field ip = — 1. This is 

Hint = {H - <?i ) H> + 2A^V 2 - a ^ty + Xip 2 ip 2 , (7.2) 
which, apart from the last term, has the form of the effective interaction for the (Reggeon) field theory of directed 



percolation |13|. With the simple rescaling 



4> = V2X/ai , <t> = VV/2A (7.3) 

one arrives at 

Hi„t = (p - ax) H + u fa 2 - 4> 2 ^j + A j> 2 4> 2 , (7.4) 
with the new effective coupling for the three-point vertices 



u 



= , M = K 2 ~ d/2 • (7.5) 



Hence u bec ome s marginal in d = 4, the upper critical dimension for DP, and indeed one would argue that the last 



term in Eq. (7.4) should then be irrelevant close to the DP transition by power counting. Thus, provided that g\ c > 0, 
we have shown that the ensuing dynamic phase transition from the inactive to the active state is characterized by the 
critical exponents of the DP universality class. Notice that both the actice and the inactive phase are characterized 
by exponential long-time behavior, as opposed to the power-law inactive phase for BARW with even to. At first 
sight countcrintuitively, the branching process accelerates the particle density decay in the inactive phase, due to the 
fluctuation-induced generation of spontaneous particle decay processes. 

Note, however, that in the bare theory the 'mass' term (the coefficient of tpij)) is negative, since fi = in that case. 
This would seem to imply that the model is immediately in its active phase as soon as o\ > 0, and hence the critical 
point at CTic = 0, which would render the rescaling ( |7.3| ) obsolete. However, this argument neglects renormalization 
effects. It may well be that the renormalized quantity — a\R is positive. As we shall argue, this is in fact the case 
for d < 2 and sufficiently small <j\ . 



23 



B. Summation of the leading singularities near two dimensions 



We are free to study the renormalization effects in either version of the theory Eq. (7.1) or Eq. (7.2). In fact, this is 
much simpler in the former, because it will turn out that the effective expansion parameter is A, and also because the 
mass term in that case is fj, + <ti, which is always positive, allowing for a perturbative expansion with a meaningful 
propagator. Even so, it is impossible to compute the renormalization to all orders. Instead, close to two dimensions we 
adopt the strategy of retaining only the most singular diagrams as e = 2 — d — ► 0, at each order in A. It is not difficult 
to see that these are given by the iterated bubble diagrams shown in Fig. ^(a,b). Notice that all the propagators are 
renormalized by a series of processes as depicted in Fig. |](c), but all these side branches may be resummed with the 
effect of replacing the mass in the internal propagator by its renormalized value. In the rest of this section, we set 
the diffusion constant D — 1, since it is not renormalized at the order we consider. The ensuing geometric sums then 
give the following simple self-consistent equations, to leading order in e, 

1 + (X/2TTe)[iJ, R + a 1R ) e / 2 

aiR = l + (\/2ne)(n R + <T 1R )-</* ' (7 - ?) 

Note that fi R + am = fi + a% in this approximation. With /i = originally, this amounts precisely to the summation 
of branched di agra ms [see Fig. Wc)] mentioned above. However, we are interested in the renormalized value of the 



DP mass (c.f. (|7J) 



(A/27re)a 1 e/2 - 1 

A R = » R - a 1R = o x \L '-^j- , 7.8 

A/27re )a t 1/2 + 1 



where we have used the above result, with /i = 0. From Eq. (7.8) we see that, for sufficiently small a\, the DP mass 
is in fact positive, indicating that the system is in fact in the inactive phase. The transition to the active phase does 
not happen until o~\ = a\ c , where, in this approximation, 

ate = (A/27re) 2 / £ . (7.9) 

Although we would not expect this result to be quantitatively correct for e = 1 , the power of A is exact and is dictated 
by dimensional analysis. Thus we expect a DP transition at a non-trivial value of o\ for all d < 2. This transition is 
driven by the fluctuations in the annihilation process. 

Since the annihilation rate is (marginally) irrelevant exactly in two dimensions, it is an interesting question whether 



it is able to drive a DP transition in this case. The single bubbles are then logarithmically divergent, and Eq. (7.8) 
is replaced by 

(A/4^)ln(A 2 /a)-l 
AR ~ (A/4*)ln(A7*) + l ' (7 ' 10) 

where A is a wave number cutoff (of the order of the inverse lattice spacing 6q 1 ). We see that indeed there is a 
non-trivial transition, at 

a lc - A 2 e- 47T/x . (7.11) 
This is one of the important results of this paper. It shows how essential it is to take into account fluctuation 



effects in studying this and similar processes. It appears to be in accord with recent detailed simulations |27|, even 
though initial work |T^| seemed to indicate that a c might be zero. Finally, for d > 2, the same approximation (which 
is, however, less justified), gives 

Ar~cti— tt — - — , (7.12) 



where 



A/(A,ai) + 1 
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In this case, as long as AA d ~ 2 -c 1, the numerator will be negative even at a\ — 0, and the system will be in the 
active phase for all u\ > 0. There does appear to be the possibility of a non-trivial transition at larger values of the 
annihilation rate, when \A d ~ 2 3> 1, but it should be noted that our approximation of summing the leading bubble 
diagrams breaks down long before then, and that it may be argued that infinite annihilation rate on the lattice already 
corresponds to \A d ~ 2 ~ 1 in the continuum theory. There is no evidence for a non-trivial transition in d = 3 from 
the simulations |17|| . 



C. RG approach 



The same problem, with m = 1, may b e at tacked using the RG loop expansion used in Sec. III. Once again, we 
work in the unshifted theory given by Eq. (7.1). The one-loop diagrams contributing to the renormalization of A, cr-y, 
and fj, are shown in Fig. ^(a,b,d), and in fact only the com bina to rial factors differ from the case m — 2 studied above. 
Defining dimensionless renormalized couplings as in Eqs. (|3.l|), (|3.2|), and in addition 



M = ZftfiK^/D 



we find the flow equations 



di{l) 
dl 

dl 
dM(l) 
dl 



= si - 
= 2M 



(l+M + Si) 2 - d / 2 

I 



(1 + M + Sl ) 2 ~ d / 2 



(1 + M + Sl ) 2 - d / 2 



(7.14) 

(7.15) 
(7.16) 
(7.17) 



Note that the effective dimcnsionless mass appearing in the denominators is now M + si. 

Adding the last two equations, we see that M(l) + Si(l) = eri(0)e 2i (equivalent to our earlier result that the 
dimensionful mass does not renormalize) . Although the remaining equations cannot be solved by quadrature, we may 
approximate their solution by setting 1 + M + s\ ~ 1 for / < Iq ~ i ln(l/<ri). Beyond that point, the couplings simply 
evolve according to their simple scaling behavior at the Gaussian fixed point. The equations may now be integrated, 
with the result for the dimensionless renormalized DP mass 



M(Zo)-s 1 (Z ) 



<j\ e 



2/„ 









^1-2 exp 


- / £(l')dl' 


) 




Jo 





(7.18) 



We may now distinguish various cases. For d < 2 and sufficiently small o\ (large Iq), the exponent in Eq. ( 7. 18] ) 
behaves like £*Iq <C 1, so that the DP mass is positive, consistent with the system being in the inactive phase. For 
large a\, however, the exponent is small and the DP mass is negative. The critical point, in this approximation, 
occurs when 



1(1') dl' = In 2 



(7.19) 



For small initial annihilation rate A we have £(V) ~ Cd\e d . This gives u\ c ~ e~ 2l ° ~ (CrfA/eln2) 2 / e , in agreement 
with Eq. (^^) apart from the factor of In 2, which may be ascribed to the crudeness of the truncation procedure 
described above. Thus, as expected, the one -loop RG equations correctly sum the most singular diagrams in 1/e in 
the weak coupling limit. However, from Eq. ( 7.19| ) we may also extract the behavior at strong coupling, because then 
1(1') ~ £* = e. In that case we find 



Clc 



-21n2/e 



(7.20) 



where, once again, the precise numerical factor is not to be taken too seriously. 



In exactly two dimensions, on the other hand, J Q ° 1(1') dl' ~ l n(l + C^XIq), c.f. Eq. (3.46), so that o\ c ~ e 



-4tt/A 

in precise agreement with the summed bubbles calculation, Eq. (7.11). The only corrections to this formula are seen 
to come from higher order loops, but, since I is no w marginally irrelevant, it should asymptotically become exact as 
A — > 0. Finally, we may use the criterion Eq. ( [7.19 ) to explore the case d > 2. In that case we may write 
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Jo £{l)dl= L wm =ln {wrm) ■ (7 - 21) 

Since now £(^o) — *• as ui — > 0, we see that there is a possibility of finding the inactive phase for sufficiently small, but 
nonzero o\ only if ln[(|e| + C*dA)/|e|] > ln2, i.e., C^A > |e|. Once again, of course, one must question the validity of 
the loop expansion in this region. However, the result does indicate that the inactive phase does not exist for nonzero 
crj at weak coupling. 



D. Larger odd values of m 



The advantage of the RG loop expansion approach outlined above for the case to = 1 is that it may readily be 
extended to larger values of to, where the identification and resummation of the most singular diagrams is somewhat 
more difficult. The only case that we shall treat explicitly here is m = 3. In that case we have a branching rate a 3 , 
corresponding to a dimensionless renormalized coupling S3, but, as before, couplings a\ and /x are generated. (Note 
that (T2 is also generated, but only at higher orders, and in any case this is irrelevant near d = 2.) Once again, the 
one-loop diagrams are simple, see Fig. |^(a,b,d) and differ from those for m = 1 only by combinatorial factors. The 
RG flow equations are 



(7.22) 
(7.23) 
(7.24) 
(7.25) 
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dl 


(1- 
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+ M + si +s 3 ) 2 " 
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6ts 3 



Notice that the mass term is now M + si + S3- As before, we see that this combination obeys a simple equation, and 
that it grows exponentially like e 21 . Thus, once again, we may make the same approximation and ignore it in the 
denominators up to Iq ~ \ ln(l/<73). The remaining equations may now be solved by quadrature. We give the results 
only for the case d — 2. If we shift the fields we now find that the DP mass is A = fi — o\ — 3<r 3 . The dimensionless 
renormalized version of this quantity at scale Iq then turns out to be 



M{l )- Sl {l )-3s 3 (l ) = a 3 e 2 



1 



12/5 



8/5 



1 + AZ /2tt (l + A? /2vr) e 



(7.26) 



Once again, this is positive for large lo (small 03), indicating the existence of an inactive phase, while it is clearly neg- 
ative for sufficiently small large a 3 . The zero must be found numerically, and occurs at Xlo/2n « 1.42, corresponding 
to 



0"3c 



,-5.68rr/A 



(7.27) 



This is expected to be asymptotically exact (apart from a possible prefactor) as A — > 0. Note that, compared with 
the case m = 1, the transition occurs at a much smaller value of the branching rate. This is because the generation of 
the process A — > now happens only at higher order in the annihilation rate A, which is asymptotically small. This 
pattern of the decrease of a rnc with increasing to continues, and is consistent with what is observed in simulations 
These calculations for larger m may also be extended to the cases d ^ 2, with similar results as for to = 1. 



VIII. BARW WITH THREE-PARTICLE ANNIHILATION 



We finally briefly discuss the case of BARW with k = 3, i.e., three-particle annihilation processes instead of the 
two-particle annihilation A + A — > which has concerned us in the bulk of this paper. In physical dimensions 
d > 1, BARW with any higher value of k will either be satisfactorily described by the mean-field rate equation (L~j). 



according to the fact that the annihilation vertex Xk has upper critical dimension d c — 2/(k — 1), c.f. Eq. (1.15), or 
simply generate the processes with lower k, see below. 
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We may of course largely follow the arguments given in the preceding sections. First, via a combination of branching 
and annihilation the branching processes with to offspring particles (rate o~ m ) now generate all the lower branching 
processes oc u m -3, ff m -6, etc., see Fig. |l^(a). Clearly, the reactions with lowest offspring number will be the most 
relevant ones, and therefore we shall essentially have to discuss only the cases m = 1, m = 2, and to = 3. For 
the latter, and equally for any m — 31, I — 1,2,..., the branching and annihilation reactions conserve the particle 
number locally modulo 3, and this symmetry precludes a fluctuation-induced generation of spontaneous particle decay 
processes, which might lead to a phase transition of the DP universality class. In analogy with the situation for even 
to in the case k — 2, there might potentially arise a new non-trivial universality class below some novel critical 
dimension d' c < 1; yet in one dimension one merely expects logarithmic corrections to the mean-field behavior. 

This leaves us with the cases to = 2 and m — 1 (or, more generally, to — 2 mod 3 and to = I mod 3). Here, the 
above special symmetry is not present, and via the processes depicted in Fig. |T^(b) and (c), respectively, the single- 
particle decay A — > is generated. For m = 2, this already happens to lowest order in the branching and annihilation 
vertices, 0(0^3), see Fig. |l^(b). In addition, the processes A — > A + A and A + A — > also become generated, 
but require diagrams of Oya^X^) and O[o\\\), respectively. This leads to all the terms required for Reggeon field 
theory to represent the effective action for the ensuing dynamic phase transition, provided U2c > (compare Sec. VII ). 
Yet as already the lowest-order diagram in Fig. [w|(b) is (logarithmically) divergent in one dimension, following the 
arguments given for the odd-offspring BARW with k = 2 we expect such a non-trivial transition to occur in d < 1, 
with the exponents of the DP universality class. In fact, as the BARW with k = 2 and m = 1 becomes generated, 
albeit to some high order in the original couplings, the DP transition may actually persist for any d < 2. For the 
case m — 1, the lowest-order diagram, 0(crfX 3 ), diverges already for d < 2 dimensions, and as again the process 
A + A — > is generated via diagrams of order 0(a\X\), this maps onto BARW with k — 2 and odd to and hence one 
has to expect a non-trivial phase transition of the DP universality class for d < 2 dimensions. 



IX. SUMMARY AND CONCLUSIONS 



We have studied branching and annihilating random walks defined by the combined processes 2 A — ► (k = 2), 
A —> (m + 1)A with both even and odd offspring number m in low dimensions, where fluctuation effects and particle 
anticorrelations are crucial. We have employed a field-theoretic representation of the corresponding diffusion-limited 
reaction problem, derived from the classical master equation. Apart from the continuum limit, the derivation of the 
action requires no further approximations (Sec. |n|). 

While the mean-field rate equation predicts the existence of an active phase only for any nonzero value o~ m of the 
branching rate (Sec. |), in low dimensions fluctuations lead to the emergence of a non-trivial inactive phase as well, 
and a dynamic phase transition at a critical value o~ c (for fixed annihilation rate A and diffusion constant D). As a 
consequence of the local parity conservation in the case of even to, which is reflected in a discrete symmetry in the 
action, both the inactive phase and the critical properties at the nonequilibrium phase transition differ drastically for 
even and odd to, respectively, as was already known from numerical simulations pll,0,ff9lp2 26 271]. Using diagram- 



matic summations and renormalization group methods, we believe we have provided at least a qualitatively satisfying 
understanding of these different universality classes || . 

For odd m and d < 2, including the borderline dimension d c = 2, fluctuations generate the spontaneous decay 
A — > 0, as well as branching processes with all odd offspring numbers smaller than m. This single-particle decay 
provides a very effective mechanism counteracting the particle production through branching, and induces a non- 
trivial dynamic phase transition characterized by the critical exponents of the directed percolation universality class, 
as has been established in simulations |I(],|lj],^6|,^7]] . In the emerging inactive phase, the spontaneous decay processes 
dominate, and thus the long-time decay is exponential, i.e. faster than in the pure annihilation model without the 
branching processes. For dimensions d > 2, we find o~ c = 0, and the system is essentially described by mean-field 
theory. Thus the values of the critical exponents display a discontinuous jump at d c = 2 (for the upper critical 
dimension of directed percolation is d c = 4). Precisely in two dimensions, we still predict a non-trivial transition with 
DP exponents, and first indications for this have recently been found (^7). Indee d, ev en the tendency that a c should 
decrease with increasing m is at least qualitatively explained by our results (Sec. VII ). 

While the case of odd m is thus well understood now, there remain a number of open questions for even offspring 
numbers. Similarly to the situation for odd to, the even to universality class is determined by the lowest possible value 
to = 2. Our analysis shows that near two dimensions, the branching rate always remains a relevant quantity (i.e., an 
active phase governed by exponential correlations), and thus necessarily a c = 0. In Sec. |l|, we have evaluated the 
critical exponents to second order in e = 2 — d, and have discussed the logarithmic corrections arising at the critical 
dimension d c — 2. On the other hand, a direct and exact computation in one dimension (Sec. IV) establishes that 
the branching rate is actually irrelevant in the limit A — > 00, which implies the existence of an inactive phase only 
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in this special situation, governed by the power laws of the pure annihilation model (see also Ref. jL7|). Therefore a 
new critical dimension 1 < d' c < 2 must exist, below which there appears a dynamic phase transition to a power-law 
inactive phase, see Fig. f?|. 

Within the framework of a truncated loop expansion to one-loop order (Sec. 0), we may describe the above 
scenario as a function of dimension d at least qualitatively; in this approximation, we find d' c = 4/3. Although the 
structure of the phase diagram seems to be correct, the ensuing one-loop critical exponents in this rather uncontrolled 
approximation are well off the numerically found values [|l7 19 2^j23|. Yet unfortunately, we were unable to improve 
on our approximations, e.g. by extending the calculations to the two-loop level. Also, the emergence of dangerous 
irrelevant variables precludes the firm establishment of scaling relations (22) relating exponents in the active phase, 
say, and those at the critical point. Although so far there has been no evidence for this, the possibility of violations of 
such scaling laws should be carefully and thoroughly reinvestigated in numerical simulations. It would also certainly 
be highly desirable to possibly design a different renormalization scheme to address the even m universality class. 
The desirability of going beyond the one-loop approximation is even clearer in the inactive phase, where the density 
is expected to decay as A/(Dnt) d / 2 , for d < 2, and Dr is a renormalized diffusion constant. To one-loop order, no 
such renormalization takes place. It is therefore not possible to establish postulated scaling laws relating the manner 
in which Dr is supposed to behave near the critical point to the other exponents defined in the active phase. The 
inactive phase is also interesting in that it is expected to exhibit a kind of spontaneous symmetry breaking as a 
result of the parity invariance of the dynamics. That is, if we start with a state with a even number of particles 
(equivalent to imposing periodic boundary conditions on the corresponding Ising model) , the asymptotic state is the 
trivial absorbing state (fully ordered in the Ising language). But if we begin with an odd number of particles this 
cannot happen: instead we expect the system to evolve to a state with a finite average total number of particles. 
Once again, one may postulate scaling laws relating the power law divergence of this number as the critical point is 
approached to other exponents, and it would be very useful to have a more systematic derivation of these, which is 
so far lacking in our analysis. The study of the nature of the inactive phase seems to be particularly interesting, as 
numerical studies appear to indicate the existence of unexpected simplifications Q . 

As discussed in Sec. (vj, an obvious generalization to N > 2 particle 'flavors' leads to qualitatively new behavior, 
governed by the exactly computable exponents of the limit N — > oo. To our knowledge, this case of BA RW w ith even 
m and N > 2 particle types has not been studied in simulations, nor have BARW with k — 3 (see Sec. VIII). 
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FIG. 1. Elements of the perturbation theory for BARW with k — 2: (a) propagator (— iw + Dq 2 + a m ) 1 ; (b,c) annihilation 
vertices A and —A; (d) branching vertex a m ; and (e) pair production vertex r. 

FIG. 2. One-loop diagrams for the renormalizations of (a) the annihilation vertex A, (b) the branching vertex a m , and (c) 
the pair production vertex r; (d) the generation of the branching process oc <r m _2 through a combination of branching cx a m 
and two-particle annihilation. 

FIG. 3. UV-singular diagrams (near d c — 2) to two-loop order for the vertex functions (a) (b) T^^,, and (c) T^. 

FIG. 4. An important set of diagrams contributing to P m +i(t), i.e., to the renormalization of the branching process, for 
m = 2. 

FIG. 5. Complete set of tree and one-loop diagrams for BARW with m = 2, for the two- and four-point vertex functions 
(a) r^, (b) T^, (c) (to first order in r ), (d) T^^, and (e) V^^. 

FIG. 6. Flow diagram s 2_d//2 vs. I for BARW with m = 2 in d = 1. For the flows depicted here, the initial values for 
the couplings are £(0) = 0.01, with (clockwise) s(0) = 0.01 (long-dashed), s(0) = 0.0005 (dotted), s(0) = 0.00025 (dashed), 
s(0) = 0.0002 (solid), s(0) = 0.00015 (dotted), s(0) = 0.0001 (dashed), and s(0) = 0.00005 (solid). 

FIG. 7. Phase diagram for BARW with even m as function of dimension. 

FIG. 8. Bethe-Salpeter equations for (a) the annihilation rate A, (b) the branching rate a', and (c) the pair production rate 
r of the N species model in the limit N — > oo. 
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FIG. 9. Most singular (bubble) diagrams for the renormalization of (a) the particle decay rate n, (b) the branching rate cti, 
and (c) the propagator (mass) /i + a for BARW with m = 1 near d c — 2 dimensions. 



FIG. 10. Diagrams for BARW with 3A — ► for the generation (a) of the branching prOC6SS OC Cm — 3 from <7 m and A3, (b) of 
the particle decay A — > for the case m = 2, and (c) for m = 1. 
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FIG. @ (J.L. Cardy and U.C. Tiiuber) 
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FIG. | (J.L. Cardy and U.C. Tiiuber) 
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FIG. @ (J.L. Cardy and U.C. Tiiuber) 
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